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We give a complete and detailed description of the computation of black hole en- 
tropy in loop quantum gravity by employing the most recently introduced number- 
theoretic and combinatorial methods. The use of these techniques allows us to per- 
form a detailed analysis of the precise structure of the entropy spectrum for small 
black holes, showing some relevant features that were not discernible in previous 
computations. The ability to manipulate and understand the spectrum up to the 
level of detail that we describe in the paper is a crucial step towards obtaining the 
behavior of entropy in the asymptotic (large horizon area) regime. 
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I. INTRODUCTION 

The identification of the microscopic degrees of freedom responsible for the black hole 
entropy is one of the key results expected from any candidate quantum gravity theory. In 
this respect loop quantum gravity (LQG) |[1]-H can claim a reasonable success because 
it accounts for the black holes degrees of freedom in a beautiful and mathematically 
appealing way |]4l [5|. Black hole degrees of freedom are described in this framework 
as Chern-Simons (CS) states residing on the "surface" of a black hole modeled as an 
isolated horizon ||6l. This leads to an interesting interplay among quantum geometry, 
Chem-Simons theory, and statistical mechanics. When the states describing a black hole 
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are counted and this number is related to the black hole horizon area a, the resulting en- 
tropy is found to be proportional to a. By taking advantage of the free Immirzi parameter 
y that labels inequivalent quantum sectors of general relativity in the connection formu- 
lation, it is actually possible to make the proportionality factor between the entropy and 
the area equal to 1/4, thus reproducing the Bekenstein-Hawking area law. The choice 
of y, initially made in the study of spherical non-rotating black holes, is universal in the 
sense that it works also for other types of black holes (i.e. rotating or coupled to some 
matter fields). As of today, there are no exceptions to this universality (see, however, 
HZl for some interesting suggestions about the role of y). The available framework also 
predicts the first subleading corrections to the Bekenstein-Hawking law which turn out 
to be logarithmic with area. These corrections are generically independent of y and in 
qualitative agreement with the ones obtained by using completely different approaches 

In addition to the successful derivation of the entropy-area law, there have been some 
recent results concerning small black holes |!9]-iTT|. In particular, a persistent "periodic- 
ity" has been found in the so called black hole degeneracy spectrum when plotted as a 
function of the areaj^ The most degenerate quantum configurations accumulate around 
certain evenly spaced values of area, with a much lower degeneracy in region between 
those values. This produces an effectively equidistant area spectrum, despite the fact that 
the area spectrum in LQG is not equidistant. This phenomenon makes contact in a non 
trivial way with the evenly spaced black hole horizon area spectrum predicted by Beken- 
stein and Mukhanov in |[12l [T3B under quite general conditions. The periodicity in the 
degeneracy spectrum leads to a striking staircase behavior when the entropy is plotted 
as a function of the area. 

In a recent paper [14] we have proposed a novel way to understand the black hole 
degeneracy spectrum by relying on number-theoretic and combinatorial methods. The 
main goal of that paper was to confirm the original results obtained in [|9l-[TT|| and extend 
them by using an improved computer algorithm based on the new mathematical under- 
standing of the problem. Once this goal has been achieved, the real challenge is trying to 
see if macroscopic black holes display the same periodicity in the entropy as the micro- 
scopic ones. To this end, one has to find appropriate exact expressions for the entropy as 
a function of the area, which are suitable to derive its asymptotic behavior. A key step 
in this direction is to codify the solution to the combinatorial problem of computing the 
number of relevant microstates as the coefficients of a formal power series expansion of 
the so called generating function (see [ jlSl [T6| for a pioneering suggestion in this direc- 
tion). In many cases this formal expansion converges in non empty open disks in the 
complex plane and admits an analytic extension allowing the coefficients to be written as 
contour integrals whose asymptotic behavior can then be studied by more or less stan- 
dard methods. The relevant generating functions for our problem were obtained and 
used to derive the expression of the entropy as a function of the area in [[T7l[T8| . 



The study of the asymptotics of black hole entropy has been attempted before |I19| . 
In this remarkable paper, Meissner gave a clever method to sidestep many of the fine- 
grained issues that we have looked at with our number- theoretic methods. There he 
obtained a closed expression for the entropy as a function of the area in the form of an 
inverse Laplace (actually Laplace-Fourier) transform. This poses some interesting ques- 

^ Here and in the following, periodicity actually means a modulation with a regular period of some growing 
magnitude. 
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tions about the compatibility of his approach and ours: Are they compatible? Is it really 
necessary to resort to the number theoretical methods of [14^ llZll? These questions have 
been answered in [18] (and the present work offers more evidence supporting the need 
to use our methods to satisfactorily understand several important issues related to black 
hole entropy). A first relevant point to highlight is the difficulty in confirming the struc- 
ture seen in [lOJ by directly computing Meissner's integrals due to the subtly oscillatory 
nature of their integrands. To this end, the number theoretical and combinatorial meth- 



ods of [|14|| are clearly better suited. The issue is then to show that they provide the same 



answer as Meissner's approach. This has been done in [T8| where the integral expressions 



of Ifl9l have been obtained from the generating functions introduced in |fT7|. Another im- 
portant issue that has been considered in [18] is the analytic structure of the integrand 
of the inverse Laplace transform that gives the entropy. The presence or absence of the 
sought for periodic structure for macroscopic black holes depends on the location of the 
poles of the integrand. The results presented in [T8| do not prove that the periodicity ex- 
ists for large black holes, but certainly do not exclude this possibility owing to the subtle 
behavior of the real parts of the poles. An alternative way to understand the periodicity 
of the entropy (and justify why it should be present) has been given in [20]. We want 
to mention here that some non-trivial "intermediate" regimes may be significant in the 
behavior of the entropy because it is not inconceivable that the sizes of real astrophysical 
black holes are actually outside the large-area asymptotic regime. This would mean that 
a non-trivial behavior could be displayed by these objects even if the asymptotics of the 
entropy is more or less trivial (i.e. just linear). Finally, the number-theoretic methods 
that we propose lead to a clear identification of the role played by the different sources 
of degeneracy in the black hole spectrum and how they explain the observed non-trivial 
behavior of black hole entropy in LQG. 

The purposes of this paper are manifold. First of all, we give a detailed and expanded 
derivation of the number-theoretic methods of [14] and show how they can be effectively 
used to understand in full detail many issues related to the behavior of black hole en- 
tropy. To avoid confusion, we consider separately the computation of the entropy in the 
Ashtekar-Baez-Corichi-Krasnov (ABCK) framework [|5l, as reinterpreted by Domagala 
and Lewandowski (DL) [|21|, and other countings proposed in the literature [ 22tl23ll . We 
also consider the novel proposal of ["24] where the authors study black hole entropy from 
first principles without using any internal gauge fixing of the SU(2) symmetry on the 
horizon. The reason to consider these other points of view is to illustrate the flexibility of 
our approach. Our second goal is to derive and discuss the black hole generating func- 
tions that codify all the information about black hole entropy as a function of the horizon 
area. We will show how they can be immediately applied to derive the asymptotic be- 
havior of the black hole degeneracy spectrum for some important subsets of the area 
spectrum (a discussion of how they can be used to get an integral representation for the 
black hole entropy appears in |[T8l ). The third goal of the paper is to provide a detailed 
account of how the periodic pattern in the black hole degeneracy spectrum arises. We 
will take advantage of the number-theoretical methods described in the beginning of the 
paper to study the interplay between the different contributions of the different combina- 
torial factors to the entropy and also to show what the most relevant area configurations 
are. We will also provide a simple heuristic argument that complements the one given 
in [|2Qj]. Another important result in this respect will be to show that the interesting be- 
havior of the black hole entropy is present even without taking into account the so-called 
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"projection constraint" (though in this case the logarithmic correction to the entropy is 
absent). Finally we give a number-theoretic prescription to label and classify the sub- 
structures ("bands") that appear in the black hole degeneracy spectrum. Summarizing, 
the main new results presented in the paper are: 

• We provide several independent methods to count the configurations selected by 
the projection constraint, in particular, we explain how the counting can be made 
by using auxiliary quantum Hamiltonians for spin systems. We also give a group- 
theoretic treatment that can be used for all the different countings discussed in the 
literature. 

• We explain how to obtain the asymptotic behavior of the black hole degeneracy 
spectrum for subfamilies of area eigenvalues characterized by a set of square free 
integers. We use then this method to show how the band structure of the degener- 
acy spectrum originates. 

• We discuss the relative importance of the different sources of degeneracy. In partic- 
ular we disentangle the contributions coming from the projection constraint from 
the remaining ones. 

• We provide a new argument that explains how the configurations contributing to 
the bands in the degeneracy spectrum can be labeled by a simple function P of 
the spin labels at the punctures on the black hole horizon. We do this by using a 
straightforward continuum approximation. 

• We introduce a new type of generating function, defined with the help of P, that 
allows us to isolate the configurations contributing to specific bands in the degen- 
eracy spectrum. 

• We give the full details regarding the degeneracy spectrum and entropy for the 
black holes with the smallest possible areas. 

The paper is organized as follows. After this introduction we devote the next section 
(section to a quick discussion of how black holes are modeled with the help of isolated 
horizons. The aim of this section is to provide a frame to assess and compare the differ- 
ent proposals, in particular those suggesting to employ a SU(2) Chern-Simons model to 



describe the quantum degrees of freedom on the horizon. Section III describes in detail 
the methods used to count the number of states relevant for the computation of the en- 
tropy the ABCK prescription according to DL. In this section we introduce several precise 
definitions (black hole configurations, degeneracy spectrum, etc.) and solve the relevant 
number-theoretic and combinatorial problems step by step. Section |W] will discuss the 
application of our methods to other countings and proposals, in particular, there is a sub- 
section devoted to the Ghosh-Mitra counting [22\ and another to the SU(2)-black entropy 
proposal of Engle, Noui and Perez | |24| . In the process we will describe several ways to 
solve the so called projection constraint. One of them is interesting because it offers some 
tantalizing hints about the possible connection of the problem of computing black hole 
entropy with conformal field theory techniques [25]. Another one makes use of a simple 
type of generating function (a fact that originally suggested that it might be possible to 
find generating functions for the exact black hole entropy (13 )• A very short summary of 
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the different schemes is presented in section |Vj Section |VI] provides a detailed account of 
the features of the black hole degeneracy spectrum, how it arises, what the most relevant 
area eigenvalues are, and the mathematical classification of the characteristic bands that 
it shows. We give here an argument, alternative to the one given in [20], that explains 
how the periodicity of the black hole degeneracy spectrum arises and we also provide 
a generating function that selects the configurations that define the bands. In particular, 
the analysis presented in this section contains new important information about the de- 
tailed features of the entropy spectrum, and introduces new techniques for the study of 
its periodic structures. All this information is valuable for the purpose of obtaining the 
asymptotic behavior of the entropy in the limit of very large (astrophysical) black holes, 
which is the main open problem in the number theoretic approach that we are follow- 
ing. Section VII is devoted to our conclusions, comments, and suggests directions for 
future work. We end the paper with several appendices. The first one (Appendix |A]) pro- 
vides a pedagogical derivation of the generating functions used in the paper. Appendix 
|B] gives a complete and explicit computation of the black hole degeneracy spectrum and 
the entropy for all the small black holes with areas smaller than 1 8 x Anyl^. (Throughout 
the paper we will denote by £p the Planck length.) Finally, Appendix |C] gives a unified, 
group-theoretic treatment to solve the so called projection constraint in any of the count- 
ing proposals discussed in the paper. 



II. ISOLATED HORIZONS AND BLACK HOLE ENTROPY 

The purpose of this section is to give a short history of the study of black holes in 
the loop quantum gravity framework and, in particular, of the entropy computations 
leading to the Bekenstein-Hawking law. Our main goal here is to provide the reader with 
the background information necessary to establish the connections between the different 
approaches and proposals, assess their degree of rigor and state of development, and 
compare their relative merits. Several excellent reviews on the subject are available (see 
for example 101261). 

The successful derivation of the Bekenstein-Hawking law is a key stepping stone in 
the quest to arrive at a working theory of quantum gravity. This explains why the 
main contenders in this area of fundamental physics -string theory and loop quantum 
gravity- have struggled to derive this result within the respective frameworksjj In both 
approaches, the first problem that must be faced is the appropriate description of a quan- 
tum black hole or, at least, a physical approximation to it. The history of the study of 
black hole entropy in LQG can be divided in two different periods separated by the land- 
mark paper of Ashtekar, Baez, Corichi and Krasnov 111]]. The main importance of this 
work lies in the fact that it proposes the key idea of describing black holes by using iso- 
lated horizons. The reason why this is so important is because a "reduction" of general 
relativity consisting of spacetimes admitting isolated horizons as inner boundaries can 
be described within the Hamiltonian formalism. This idea is similar in spirit, but not 
equal, to the study of the quantization of the symmetry reductions of general relativity 
that gives rise to the popular mini and midisuperspace models. For both types of systems 
a subset of the possible gravitational field configurations is chosen by imposing some re- 
striction on the configuration space of metrics. In the mini and midisuperspace models 



^ An interesting discussion of the relative merits and problems of both approaches can be found in ETI . 
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this is accomplished by restricting the allowed configurations to metrics satisfying some 
kind of symmetry principle whereas in the case of black holes the requirement is that the 
allowed metrics must have an isolated horizon. 

During the period preceding the appearance of H, several authors made interesting 
suggestions that have played a key role in the development of the standard ABCK model. 
In fact, it is fair to say that many of the ideas that have been instrumental in the currently 
accepted framework made their appearance during those years. We want to mention 
explicitly the work of Smolin in |l28i, where he discusses, among other things, the impor- 
tance of considering the quantization of the gravitational field in spacetimes with inner 
boundaries of finite area and the important role of the Chern-Simons theory and quan- 
tum groups. After this paper, it is mandatory to highlight the works of Krasnov 11291 130II 
and Rovelli 13111321 . The first author gives several important insights related, in particu- 
lar, to the issue of the distinguishable character of the spin network labels on the horizon, 
the types of spin networks that are relevant to describe black holes and the concrete role 
of the CS theory to describe horizon degrees of freedom]^ Rovelli, on his part, made 
a direct counting and used simple combinatorial methods to arrive at an approximate 
Bekenstein-Hawking law (in the sense that the linear growth of the entropy as a function 
of the area was found but the proportionality coefficient was not fixed). He also insisted 
on the role played by the distinguishability of punctures on the black hole horizon. 

The currently accepted treatment of the problem is given in [|4n5J (see also the review 
by Corichi [26J). The starting point is the Hamiltonian description of spacetimes with 
isolated horizons in terms of Ashtekar variables. A crucial issue is the treatment of the 
inner boundary and the consequences that its introduction has for the final Hamiltonian 
description of the model. The following facts are particularly relevant. 

• The isolated horizon condition translates itself into a matching boundary condition 
involving the puUback of the Ashtekar connection and the tetrads to the spherical 
slices of the inner boundary. 

• The inner horizon boundary conditions force the SU(2) connection A^, to be re- 
ducible on the boundary [4J (here the indices i and j are su(2) indices). This means 
that there exists an internal vector field r\ defined on the spherical section of the 
isolated horizon, satisfying the condition dAt^ = dr^ + A"^ , r' =0. A choice of this 
internal vector [that amounts to a partial gauge fixing of the SU(2) symmetry] is 
used to explicitly take into account the reducibility of the connection on the bound- 
ary. The U(l ) invariance on the isolated horizon is guaranteed by the projection 
of the matching conditions. 

• The symplectic structure consists on a volume part (the same as in the usual models 
without boundary) and a boundary part corresponding to a U(l )-CS theory. This 
fact strongly suggests that the Hilbert space appropriate to describe the system is 
a tensor product of a U(1 ) Chem-Simons Hilbert space and a LQG volume Hilbert 
space. 

• In the quantum formalism the r^-projected matching condition is quantized and 
used to select the physical quantum states. The solutions to this condition automat- 
ically satisfy the Gauss law. This means that it is treated as a first class constraint as 
far as quantization is concerned. 



^ It is interesting to mention at this point that the original proposal made use of a SU(2)-CS model. 
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The reducibility condition of the connections on the horizon is a somehow unexpected 
feature of the present scheme. Several authors have explored the physical consequences 
of heuristically forgetting about this fact and retain a SU(2)-CS model on the horizon 
l|23l |33| , whereas other authors have addressed the development of a SU(2) invariant 
formalism from first principles Il24l l34ll . The possibility of using a SU(2) description is 
an interesting topic. In fact there is a recently published but old proposal by Krasnov 
and Rovelli [35J to describe quantum black holes -without the use of any semiclassical 
approximations- that lead to results very similar to the ones presented in [24, 34 J. 

Within all these models it is possible to obtain the Bekenstein-Hawking law by appro- 
priately adjusting y. There are also logarithmic corrections to the entropy that play an 
important role because they are specific predictions of the different proposals that can be 
used, in principle, to choose among them. The combinatorial problems associated with 
the computation of the entropy for all these schemes are rather similar in nature and can 
be approached with the techniques that we describe in this paper. In that respect, our 
methods are quite robust and well adapted to the nature of the combinatorial problems 
that appear, so we expect them to play a significant role also in the understanding of 
dynamical aspects of black hole physics such as Hawking radiation or black hole evapo- 
ration. 

III. THE ABCK QUANTUM ISOLATED HORIZON AND THE DL COUNTING 

The study of black holes within LQG f5\ makes use of the isolated horizon concept 
(see |'6l and references therein). In this framework the horizon is introduced as an in- 
ner boundary of the classical spacetime manifold. Several conditions are imposed on it 
to guarantee that the relevant physical features of a black hole are captured, in particu- 
lar its thermodynamical behavior. Spacetimes of this class admit 3-dimensional, spatial, 
partial Cauchy surfaces. Each of them is bounded by a topological 2-sphere, that we will 
refer to as the horizon. A Hamiltonian description for this sector of GR is available and is 
the starting point for canonical quantization [36] . The details of the canonical description 
strongly suggest that the appropriate Hilbert space should be built as a tensor product of 
bulk and horizon Hilbert spaces, ^^^^ = ^Hor ® 3^bu1- In this approach, it is natural to 
adapt the bulk Hilbert space JCbui to the presence of an inner boundary. As usual in LQG, 
it is convenient to use a spin network basi^ A particular bulk spin network can pierce 
the horizon or not. If it does, it will do so at a finite number of points that we will refer to 
as punctures. These punctures have a distinguishable character, as a consequence of the 
action of general diffeomorphisms over the horizon surface [5J. The punctures carry two 
quantum numbers (j, ra). The first one is just the spin j-label of the edge that defines it. 
These spin quantum numbers allow us, in particular, to calculate the horizon area accord- 
ing to the standard prescription of quantum geometry. The other quantum numbers are 
spin components, m, defined with the help of the preferred 5u(2)-internal vector field r*^ 
on the horizon (see section |ll|). These labels play an important role in the implementation 
of the quantum boundary conditions. Explicitly, the bulk Hilbert space CKbui is spanned 
by states of the form 

1(0), •• • )bu1, . . . , Ktui, ji, . . . , raN, Jn), • • • )bu1 , • • • (ULl) 

^ A spin network is an oriented graph embedded in a spatial section of the spacetime, carrying a SU(2) 
irreducible representation labeled by a spin number j on each of its edges, and a gauge invariant operator 
(intertwiner) linking incoming and outgoing representations at each vertex [3|. 
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where the half integers ji G \N label irreducible representations of SU(2) and 

mi G{-jh + (111.2) 

are spin components in the direction of r\ The numbers mj and ji represent the quantum 
degrees of freedom of the bulk geometry "close" to the horizon, and the "■ • • " in the bulk 
state |(mi , ji , . . . , raN, Jn), ■ ■ ■ )bu1 refer to bulk degrees of freedom away from the horizon. 
Finally 1(0), ■ ■ ■ )bu1 denotes the states corresponding to spin networks that do not pierce 
the horizon. 

The surface horizon Hilbert space CH^Hor described by a U(1) Chern-Simons theory. 
The level k G N of this quantum CS-theory gives rise to a prequantized value = Anyly k 
for the area of the isolated horizon. In 3^Hor' the U(1)-CS basis states over the punc- 
tured sphere |(ci, . . . , CN))Hor ^^"^ characterized by arbitrarily long (ordered) sequences 
(ci)[ii — (ci, . . . , On) of non-zero congruence classes of integers modulo k. Each of the 
Ci can be thought of as an integer number belonging to {1 , 2, . . . , k — 1} that labels the 
quantized deficit angle 47tci/ k associated with the I-th puncture. If the isolated horizon 
has a spherical topology, there is an additional restriction over the total curvature, which 
translates into a condition 

^ Ci = (mod k), 
I 

for the congruence classes appearing in a given sequence. This condition can be inter- 
preted as the quantum equivalent of the Gauss-Bonnet theorem. These Ci labels turn out 
to be related, through the quantized isolated horizon boundary conditions, to the label 
mi in the ji representation of the spin network edge, | (ji, mij^i^ , ■ ■ ■ )bu1 G 3^bu1/ piercing at 
the corresponding puncture, according to 

Cj = — 2mi (mod k) , 

with mi G {— ji, — ji + 1 , . . . , ji}. This restriction on the form of the basis states |(ci))^or ® 
Kji, mi), ■ ■ ■ )bu1 of CK^jj, is the quantum counterpart of the isolated horizon boundary con- 
dition. 

The area induced on the horizon by a given bulk state is the eigenvalue of the standard 
area operator 

a^«'=(ji, mi) = ?>nyil Y_ VJi(Ji + 1) , (ni.3) 

I 

given in terms of the jplabels of the punctures. For a fixed value of the black hole area, 
there are different combinations of horizon labels compatible with it (i.e. the area eigen- 
values are degenerate). The quantum states of the horizon belonging to ^K^^j. and "com- 



patible" with the area (III.3) are precisely the ones responsible for the black hole entropy. 



Some important comments, relevant to the definition of the entropy, are in order now: 

• Compatibility of areas: Two different types of areas associated with the horizon 
are relevant here. The prequantized area Qk necessary for the quantization of the 



CS theory and (III. 3). The does not belong to the spectrum of the area operator, 
i.e. 

Qk = AnyllK ^ a^'2^(ji, mi) = Myll }^ Vji(ji + 1). 
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This forces us to introduce a suitable, and physically sensible, notion of compatibil- 

• Entropy: In Quantum Statistical Mechanics the standard definition of the entropy 
in the microcanonical ensemble makes use of an energy interval 1381 . In vacuum 
general relativity, the non-rotating and neutral black holes are fully characterized 
by their areas. Therefore, it is convenient to use this geometrical feature in the dis- 
cussion of their statistical mechanics. This is explained at length in ||29ll39]l . Here, 
then, we introduce an area interval [a^ — 5, Qk + 5] in order to define an appropri- 
ate statistical ensemble. This helps to solve the problem mentioned in the previous 
item because, once the prequantized value of the area is fixed, there are always bulk 
quantum states with area eigenvalues belonging to [a^ — 6, + 6] for reasonable 
choices of 6. The entropy is obtained by tracing out the bulk degrees of freedom 
to get a density matrix that describes a maximal entropy mixture of surface states 
with area eigenvalues in the previous interval. The value of the entropy is com- 
puted by counting the number of allowed sequences (ci) of non-zero elements of 
Zk satisfying Ci + ■ ■ ■ + Cn =0, such that Ci = — 2mi (mod k) for some permissi- 
ble spin components [mi]. Here permissible means that there exists a sequence of 
non-vanishing spins (ji) such that each mi is a spin component of ji and 

aK - 6 < a^'2"(ji, mi) = Syry^i Y. + < aK + 6. (III.4) 

I 

The counting of horizon c-labels is equivalent to the determination of the dimension 
of the Hilbert subspace of ^K^oj. that represents the black hole degrees of freedom. 

• Thermodynamical limit: In the thermodynamical limit, for standard statistical me- 
chanical systems, the choice of energy interval introduced in the definition of the 
microcanonical ensemble is irrelevant II38II , i.e. it is equivalent to consider [E — A, E + 
A], or [E — A, E], or [E, E + A], or even [Eq, E], where Eq is the minimum energy of the 
system. Moreover, only in the thermodynamical limit the entropy is a smooth func- 
tion of the energy with the possible exception of points related to phase transitions. 
This smoothness is crucial to define the derivatives appearing in the definition of 
thermodynamical magnitudes. In the present case, one should follow a similar path 
(with the area playing the role of the energy). 

The possibility of taking an appropriate area interval has been used to simplify the 
actual computation of the entropy. By using an interval of the form [0, a^], Domagala and 
Lewandowski [21J proved that the black hole entropy can be obtained according to the 
following prescription involving only the raj bulk labels: 

The entropy S<(aK) of a quantum horizon of classical area Uk, according to Quantum Geom- 
etry and the Ashtekar-Baez-Corichi-Krasnov framework, is 

S<(aJ=log(1+N<(aJ), 

^ In this respect we want to point out that there is another possible choice |l37| for the horizon area operator 
in LQG corresponding to an evenly spaced area spectrum. In this case the Bekenstein-Mukhanov IIT2l[T3l 
scheme is realized from the start and there is no need to introduce an area interval at this stage. 
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whereN<{ati) isthenumber of all the finite, arbitrarily long, sequences (mi, . . . , trn) ofnon-zero 
half integers, such that the following equality and inequality are satisfied: 

N N 

^mi = 0, ^VlTni|(|mi| + l)<^. 
The extra term 1 above comes from the trivial sequence. 

Notice that, in the ABCK prescription, the entropy S< (a) is defined only for area values of 
the form a= a^. However, we will extend the definition to arbitrary values of a e [0, oo) 
by just requiring that 

|:Vlm,l(|m,| + l)<3^. 

We introduce this extension because we think that the detailed form of the area spectrum 
predicted by loop quantum gravity may play a prominent role in gravitational systems 
for which area is an important observable quantity. In the following, unless stated other- 
wise, we will write areas in units of Aixyl^. 

The rest of this section is devoted to describing methods based on number theory 
that are useful to understand the structure of the black hole degeneracy spectrum and 
entropy. The counting of the allowed sequences according to the previous definition is 
conveniently performed in four successive steps: 

Step 1. Fix a value for the area a and obtain all the possible choices for the half integers |mi| 
compatible with the area, i.e. satisfying 

N 

^V|mi|(|mi| + l) = |. 

1=1 

Notice that, at this point, we are considering the possible choices of absolute values 
of the spin components as the elements of a multiset (and hence there is no ordering 
of the labels). In other words, at this stage we only want to know how many times 
the spin component 1 /2 appears, how many 1 's appear, and so on. 

Step 2. Count the different ways in which the previous multisets can be reordered. 

Step 3. Count all the different ways of introducing signs in the sequences of positive half- 
integers I mi I obtained in the previous step in such a way that the projection con- 
straint Xjli mi = is satisfied. 

Step 4. Repeat the same procedure for all the area eigenvalues smaller than a and add the 
number of sequences obtained in each case. 



III.l. Step 1: The Pell equation 

The first step can be thought of, in fact, as a characterization of the part of the spectrum 
of the area operator relevant to the computation of black hole entropy (regarding, in 
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particular, the degeneracy of the different area eigenvalues). Let us start by introducing 
the positive integer variables kj := 2|rai| and write 

N N kmax 

2Y_ Vlmi|(|mi| + 1) = a ^ V(ki + 1)2-1 = a ^ N^y/J]^+W^ = a (III.5) 

1=1 1=1 k=l 

where the non-negative integers Ni^ (that will be allowed to be zero) in the last sum tell 
us the number of times that the label k/2 G N/2 appears in the sequence (|mi|) = (ki/2). 
Also, we denote as k^ax = k^ax ( a) the maximum value of the positive integer k compati- 
ble with the area a. At this point we are interested in finding out the multisets mentioned 



above, i.e. all the sets of pairs {(k, Ni,) : k G N, G N U {0}} such that equation ( |III.5[ ) 
is satisfied. Notice that in the description provided by the multiset, we can restrict our- 
selves to list only the values of k that do actually appear (i.e. those for which N]^ ^ 0). 
We want to point out now a simple but important fact. By using the prime factor de- 
composition of k(k + 2), it is always possible to write A/(k+ 1)^ — 1 as the product of an 
integer times the square root of a square-free positive number p^. Square-free numbers 
are integers that are not divisible by non-trivial square numbers (that we enumerate as 



pi = 2, p2 = 3, p3 = 5, and so on). Then, equation (III.5) tells us that the area eigenvalue 
a must be an integer linear combination of square roots of squarefree numbers, and have 
the form 

1-max 



i=l 



or else equation (III.5 1 cannot be satisfied. This leads then to the following condition 

l^max i-max 



k=l i=l 



As a preliminary step to solve this equation, we want to separately consider each of the 
square-free numbers pi and find out the possible values of k such that A/(k+ 1)^ — 1 is 
an integer multiple of y/pl, i.e. we first solve the equations 

in the two unknowns kt and yi (here the label 1 refers to the square-free number p^). These 
are equivalent to 

x[-pty? = 1, (IIL7) 



where we have written := k^ + 1. For each fixed square-free integer pi, (IIL7) is the 
well known Pell equation BOB , whose general solution can be found in the following 
way. Obtain first a so called fundamental solution [x],y]] with the smallest positive value 
of x by using continued fractions as explained in ||40l . Once this solution is known the 
(infinitely many) remaining ones are given by (xt,iji) := (x^^,ij^), a G N, where 
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from which we get the solutions to the original equation that we will label as (Tc^,^^). 
For instance, for the first square-free pi = 2, the previous sequence starts as 



(2,2), (16, 12), (98, 70), (576,408), . . . (see section VI.l for more details). 



Once we know the values of that can contribute when a given appears in the 



value of the area, equation (IIL6 1 can be written as 

l^max i-max CO 1-max 

Y_ ^Wi^ + 1 )^ - 1 = L L ^KyicVvi = Y. q^v^ • 

k=l i=l a=l i=l 

Now we make use of the fact that the square roots of the square-free numbers are linearly 
independent over the rationals (and hence also over the integers) to show that the previ- 
ous equation is equivalent to the following system of imax uncoupled, linear, diophantine 
equations 

oo 

Y_ y = = 1 > • • • , w. (111.8) 

a=l 

Several comments are in order now. The first is that, for a fixed value of the area a 
(necessarily an integer linear combination of ^/pl with a finite number of coefficients q^) 
only a finite number of labels k]^ come into play. Second, it may happen that some of these 
equations admit no solutions, in that case a does not belong to the spectrum of the area 
operator (this happens, for instance, for a = Vl). Finally, when they can be solved, their 
solution tell us exactly what the allowed values for k are and the number of times Ni^ that 
they appear. This construction identifies the set of allowed configurations C(a) consisting 
of all multisets c — {(kj^, N^i )} associated with a value of the area a = ^.i ^iy/Vi- 



III.2. Step 2: The reordering degeneracy (r-degeneracy) 

Up to this point we have found the number of all the possible choices for the absolute 
values of the spin components |ra| = k/2 compatible with a given value of the area a to- 
gether with their multiplicities Ni^. Each of the configurations c G C( a) can be represented 
in the following schematic form 

k k kmax kma 




> T 5 • • • > T > 



2' ' 2 
\ 



dr(c) := 't^l^::' . (111.9) 



where, in the previous representation, if a particular is zero then there are no k/2 
terms. The number of different sequences (ki/2) obtained from each configuration c by 
reordering its elements is given by the multinomial coefficient 

In the following we will refer to d.r(c) as the r-degeneracy of the configuration c. We will 
also define the r-degeneracy associated with a given value of the area o. — Y.i ^iy/Vi 

D,(a) = Y_ ^r(c) . (III.IO) 

As we will see, this quantity plays a central role in the appearance of the black hole 
entropy structure found in IITOll . 
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III.3. Step 3: Solving the projection constraint (m-degeneracy) 

Once we have identified all the possible sequences (|ttli|) = (ki/2) of positive half- 
integers satisfying the area condition, we have now to introduce signs in each rai and find 
out how many of the resulting sequences (raj) satisfy the so-called projection constraint 

N 

^ mi = . (III.ll) 

1=1 

There are several approaches to solve this problem that we will describe here. The rea- 
sons to look at these different ways to solve the projection constraint are the following. 
First, some of them are specially suited to be used in the computer algorithms that we 
have employed in the actual black hole entropy computations. Second, the solution in 
terms of generating functions is a preliminary step towards the obtention of the black 
hole generating functions that we define below. Finally, some of the solutions suggest 
intriguing connections with other interesting problems (such as conformal field theories 
as we will mention briefly in Appendix |C]). Let us consider them one by one. 

The partition problem: The problem of finding all the possible different ways to "sprin- 
kle" the signs on a fixed sequence (ki/2) has already been considered in the literature. In 
fact, it has a proper name: the partition problem, that can be stated as follows. Given a 
sequence (ki/2) = (ki/2, k2/2, . . . , kN/2) of N real numbers (positive half-integers in our 
case) find all the different partitions of {1 , 2, . . . , N} = 3\f+ U 3\f_ such that 

}^ki- }^ki = 0. 

Here we will solve the following slightly different problem: given M G Z/2, find out 
the number of different ways to partition {1 , 2, . . . , N} in such a way that the following 
condition holds 

^ ki - ^ ki = 2M . (III.12) 

ieN+ ie:M_ 

The answer to this question is known and can be found, for example, in ||4T|. It is given 
by 

L- 1 N L-1 N 

f_ Y_ e-4m«M/L -Q cos(27t£ki/L) = ^ ^-^nilM/L "Q cos(47r£|mi|/L) (III.13) 

£=0 1=1 e=o 1=1 

where L is a conveniently chosen integer. As the result is independent of this choice (as 
long as L is big enough, see below) we will take L := 1 + 2M + Y-^=^ ^i- 



The formula (III.13 1 is obtained in the following way [4T1|. For a given sequence of pos- 
itive integers (ki , k2, . . . , k^j) G N"^ and M. G Z/2, let us define the auxiliary Hamiltonian 
operator H : C^®^ ^ C^®^ as 

N 

H = 2M-_^kiCTj^', 
1=1 
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where the operator aj = 1 (g) ■ ■ • (g) 1 (g) cr^^' ® ^ ® ■ ■ ■ ® ^ acts as the 0'^' Pauli matrix on 
the I* C^-factor of C^®*^ and trivially on the others. This operator satisfies 



L-l 



L-1 



/2m£ / 



-}^Tr(^exp(-2mm/L)j =-}^ Y_ exp (^— }^ kjai - 2M 

" " f=0 {tTi=±l} ■ " 



{=0 



1=1 



When ^i^^i - 2M = we have 



1=0 1=1 



On the other hand, when kiffj — IM ^ we can write 



L-l - .„ N 

L-p(^(Z 

e=o 1=1 



kiCTi - 2M 



1 



1 - exp (inif Z.i=^ l^iCTi - 2M 



exp (27ii(^X.Z^ ^i^i - /L 



0. 



In this last equation we have used the fact that Y-^=] ^I'^i ~ ^ ^ ^^i*^ L = 1 + 2M. + 
Y_i=) to guarantee that the denominator in the previous expression never vanishes. 
This way we conclude that 



-^Tr(^exp(-2mm/L)) = 5(2M,^kiai) 

f=0 {CTI=±1} I 



(III. 14) 



i.e. the trace above counts the partitions such that (III.12I holds. The left hand side of 
equation (III.14I can be explicitly computed as 

, L-l ^ L-l N 

- Y_ Tr exp(-2m£H/L)) = - ^ g-4m«M/L "Q (exp(2m£kia(3Vl-)) 

e=o 1=0 1=1 



JIN - ■ 



1=1 

N 



,-4m«M/L 



]^cos(27r£ki/L) . 



(III. 15) 



1=1 



By comparing ( III.14[ > and (III. 15) it is clear that ( III.13[ > provides the number of solutions 
to the partition problem. It is important to notice that equation ( III.14[ > implies that ex- 
pression (III. 13) is zero whenever no solutions to the projection constraint can be found. 

Generating functions: Another method to solve the projection constraint ( III.12[ > is 
based on the use of a suitable generating function. Let us consider an ordered, finite, 
sequence (ki/2) of positive half integers of the form k/2 with multiplicities given by 
N]^ > 0, and take the following function of the variable z (a Laurent polynomial) 



Yl[z'''+z-''')=Yl[z'' + z 



(in. 16) 



1=1 



By expanding it in powers of z, it is easy to see that the coefficient of the power z^*^, 

[z''^]Yl{z^ + z-^r\ 
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is precisely the number of different ways to distribute signs among the elements of the 
multiset in such a way that the sum of all the elements k /2 of t he sequence equals M. 
In particular, if we look for the constant, i.e. [z°], term in (III. 16 1 we get the number of 
solutions for the projection constraint that can be built from a given sequence. 

A convenient way to extract this type of information is by using Cauchy's theorem. 
This allows us to extract these coefficients (for a given M) by computing the contour 
integral 



Z + Z J 



1 

2ni 



dz 

,2M+1 



(III.17) 



where C is an index-one curve surrounding the origin. By choosing for it a unit circum- 
ference, the previous integral can be written in the following useful alternative form 



k>Nv 



2N-1 

71 
2N-1 

71 



cos(2Me) ]^cos^MT^0) de 

' k 

N 

cos(2Me)]^cos(kie)de, 



(III.18) 



1=1 



where N — IM^ denotes the number of elements in the sequence (ki/2). It is important 
to point out that (III.18I and ( III.13[ > provide the same answer to the counting problem. 
The advantage of this solution is that it is specially appropriate to find generating func- 
tions for the black hole entropy as will be shown later. In Appendix |C] we provide a 
different procedure to solve the same problem based on elegant group theoretical meth- 
ods [compare equations ( III.18[ > and ( C.3[ >]. 



The definition ofm-degeneracy: Given a value of the area a = qt^/pt, we will define 
the ra-degeneracy d°^(c) of the configuration c G C(a) as 



2 



N-l 



71 



1=0 X OL 



where N = N (c) = L« ^k^ and L = L(c) = 1 + Li L« K^^.v 

Black hole degeneracy spectrum: Once we have closed expressions for the r and ra de- 
generacies, we can define the hlack hole degeneracy associated with a given value of the 
area a = Li liv^ 

~ (III. 19) 



D 



DL/ 



a 



Y_ dr(c)d°Hc). 



cee(a) 



We will refer in the following to D^^(a) as the hlack hole degeneracy spectrum for the DL 
counting. According to the previous discussion, for every value of the area a, D^'^fa) 
gives the number of sequences of non-zero half-integers satisfying the two conditions 



N 

1=1 



0, }^ V|miK|mi| + r 



1=1 



a 
1' 
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The procedure described above can be efficiently implemented in a computer to explic- 
itly obtain the black hole degeneracy in terms of the area. In order to see the fundamental 
structure of the degeneracy, we can plot the black hole degeneracy D°^(a) versus a. This 
is shown in Fig. [T} We can see that the result obtained in IITOB by a brute force analysis is 
reproduced. Specifically, we can see that the number of sequences D^^[a] is distributed 
forming a "band structure" in terms of the area. Peaks of degeneracy appear in an evenly 
spaced fashion, interspaced with regions where the degeneracy is several orders of mag- 
nitude smaller. This gives rise to an effectively equidistant spectrum [10 J. 




FIG. 1. Plot of the black hole degeneracy spectrum D^^ (in units of 10^^) in terms of the area (in 
units of 47ty£p) for a range of area values. The periodicity can be traced all the way back to the 
smaller values of the area. 



In the process of understanding the structure of the black hole degeneracy spectrum 
(III.19) shown in Fig. [ij the m and r degeneracies play different roles. In order to dis- 
entangle them it will be useful to consider an auxiliary description where the projec- 
tion constraint is ignored. This means that, once we have identified all the possible 
sequences (|mi|) = (ki/2) of positive half-integers satisfying the area condition, we in- 
troduce signs in each mi without any additional restriction. For this auxiliary problem, 
the ra-degeneracy of a configuration is given by 



i a 
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and the black hole degeneracy becomes 



(III.20) 



cee(a) 



We will sometimes refer to D^'^fa) as the "black hole degeneracy without projection con- 
straint". 



In order to take into account the inequality appearing in the DL-prescription for the 
black hole entropy, one can repeat the previous procedure for each of the relevant area 
eigenvalues a' < a and add up the resulting black hole degeneracies D^^[a']. For a hy- 
pothetical equally spaced area spectrum this task could be easily accomplished by using 
generating functions (as has been done, for example, in ||T5ll37| |). However, for the general 
case that we are considering here, one has to resort to more complicated methods based 
on the use of functional equations or a combination of generating functions and integral 
transforms flS', 19]. It should be pointed out that, as the area spectrum is a countable 
set, it is possible in principle to build the sequence (D^'^fajf'^)) consisting of the values 
of the black hole degeneracy corresponding to the n^^ area eigenvalue (ordered in such 
a way that < a^-^ ), write a generating function for it, and perform the sum by the 
same method used in [115111611 . In practice, however, this is a difficult problem due to the 
(current) lack of a simple enough algorithm to obtain the value of the n* area eigenvalue 
in terms of n. This is the reason why we are forced to use Laplace transform methods 
to complete this last step. 

In order to fix ideas, let us consider the sequence [a\^) of eigenvalues of the (relevant 
sector of the) area operator and some number sequence ((3n) related to them (for exam- 
ple the black hole degeneracies corresponding to the areas ajf^). We want to solve the 
following problem: given an area a, compute the sum 



We will take into account now the following two factq^ 

1. £(5(a-(x),s) =e-«^fora>0. 

2. If f (s) = /C(F(a), s), then the Laplace transform of F(a') da' is simply s^^f (s). 



III.4. Step 4: Adding up 




The solution when a does not belong to the area spectrum is trivially given by 



[n:a\ 




(in.2i) 



^ Here and in the following L[T{a), s) denotes the Laplace transform, expressed in the variable s, of the 
fxinction F(q). Also (f (s), a] denotes the inverse Laplace transform of the fimction f{s) in terms of 
the variable a. 
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These properties allow us to rewrite (111.21 ) in the form 



Y_ (3n = £"^ (s-^}^(3nexp(-<«'=s),a 

For the cases in which pn = D°^(a|f^) or |3n = 0^^(05^^) the sum 

^(3nexp(-<«'=s) =:P(s) 



(III.22) 



(III.23) 



neN 



can be conveniently obtained from the generating functions introduced in Appendix A 
The idea is to take first the generating function given in [17J (see Appendix A for details) 



00 00 



Z, Xi,X2, 



1 - ^ ^ f + Z IxK" 



(III.24) 



1=1 a=l 

or, in the case of ignoring the projection constraint (by putting z = 1 ), 



00 00 



-1 



G°^(Xi,X2, 



G°Hl,Xi,X2, 



1 - 2 ^ ^ xf" 



(III.25) 



i=l a=l 



and perform the substitutions xi = e^vW^ ^2 = e^^^, and so on. It is important to point 
out that, in order to deal with the projection constraint, we need to introduce an extra 
variable z. This means that, in this case, we will not directly get the function P(s) appear- 
ing in (III.23) but rather a function of P(s, z) such that the inverse Laplace transform of 
P(s, z) in the variable s is a Laurent polynomial in z, with area dependent coefficients, 
whose constant term gives the desired sum. In practice this requires the computation 
of a contour integral in z or, equivalently, an inverse Fourier transform in the additional 
variable w defined by z = e^^ [18] . By following this procedure we find 



i=l a=l 



The exponentials e ^^^v^ appearing in this function can be simplified if we use the Pell 



equations (III.7I, y'^,^ = y^KdKc + 2), to get 



pDLfs,a)) = (^1 - 2 ^ ^ g-sVKk(Kk+2) cos(cuki,))" =(^1 - 2 ^ e-^V^*^ COS cuk 

i=l a=l k=l 



-1 



Now, by using (III.22I we obtain the following expression for the entropy as an inverse 
Laplace transform (and an additional inverse Fourier transform to deal with the projec- 
tion constraint) 



expS<(a) 



1 



(27r)H 



■xo+ioo 



Xq— 100 



00 ^ 

(1 - 2 ^ e-'V^^^ COS cuk) ' e"^ ds dcu , (in.26) 

k=l 



where xq is a real number larger than the real part of all the singularities of the integrand 
in the previous expression. This expression is valid for those area values a > that do 



S<(a] 
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a 



FIG. 2. Plot of S< (a) in terms of the area (in units of Anylp). The computation of S< (a) has been 
done by using the algorithm based in number-theoretical method discussed in the paper. 



not belong to the spectrum of the LQG area operator and, in particular, it gives the exact 
value of S< (Uk) for the CS prequantized values of the area. For areas of the form the 



a 



LQG it 



previous formula gives the arithmetic mean of the left and right limits when a 

The expression ( III.26[ > can be used to study the asymptotic behavior of the entropy, 
whose exponential growth as a function of the area is explained by the presence of the 
pole in the integrand of (III.26I with the largest real part |[T9|. This pole determines, in 
particular, the value of the Immirzi parameter y = 0.237 . . . that must be chosen in order 
to reproduce the Bekenstein-Hawking law. An additional logarithmic correction of the 
form 

-llog(a/£^p) 



can also be derived from (III.26). When the projection constraint is ignored, the entropy 
is given by 



expS<(a) 



1 

27ri 



'Xfl +100 



Xq— loo 



,-1 



1 -2^e-^V' 



-s./k(k+2) 



-1 



e'^'^ ds , 



(III.27) 



k=l 



from which it is easy to see that the linear behavior of entropy with area and the value of 
y are unaltered. However, there are no logarithmic corrections in this case. 

It must be pointed out that, from a purely numerical point of view, ( III.26[ > and (III.27) 
are rather bad because of the inherent difficulties associated with numerically computing 
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improper integrals. Nonetheless, it can be seen that the formula (III.26I gives the ex- 
act result for the lowest values of the area spectrum by comparing to a direct numerical 
computation (that is impractical to extend for larger areas). Entropy computations are 
carried out in practice by using the exact values for the black hole degeneracy spectrum 
obtained by using the exact combinatorial methods described above and adding up. A 
sample result for S<(a) is shown in Fig. |2| which displays a characteristic step modula- 
tion superimposed to a linear growth. This behavior may be present for large areas as 
a consequence of the fact that the analytic structure of the integrand in ( III.26[ > is rather 
complicated and, in particular, the (real) pole that determines the growth of the entropy 
and fixes the value of y is an accumulation point for the real parts of the poles in the 
integrand (Hi. 



IV. OTHER APPROACHES 



In the preceding sections we have discussed the computation of the entropy within the 
framework developed in fillSl and expanded in ||T9ll2T| . These papers have motivated 
other authors to propose modifications in which some elements differ from the original 
approach. In particular, Ghosh and Mitra have put forward a modified counting scheme 
where the set of spin labels of horizon states is expanded to include ji and rai labels [22|. 
There is another series of papers by Kaul and Majumdar that explore the possibility of 
using a SU(2)-CS theory to describe the horizon degrees of freedom ||23ll33ll . We think that 
all these proposals have, at least, a heuristic value. Finally, a concrete way to show that 
the horizon degrees of freedom are indeed accounted for by a SU(2)-CS theory is given 
in |l24l |34|. As these schemes differ in their predictions for y and /or the logarithmic 
corrections for the entropy, it should be possible in principle to falsify some (or all) of 
them on physical grounds. Here we take a pragmatic approach motivated by our desire 
to show that all these models can be rather easily handled by the techniques discussed 
in the present paper. This shows that our methods are a sharp tool that can be used to 
derive definite predictions within a broad class of LQG inspired models. The result of 
the analysis of the different schemes can be summarized by saying that they all lead to 
the Bekenstein-Hawking law (for different choices of y), they predict similar logarithmic 
corrections (though their coefficients differ among the proposals) and all of them display 
the interesting substructure found in [,9, JOJ. 



IV.l. The Ghosh-Mitra counting 

We will start by considering the alternative proposal put forward by Ghosh and Mitra 
(GM) in [22J and showing how it can be handled by the same methods that we have used 
for the ABCK approach. We will then analyze the results and compare both models. It 
is important to highlight at this point that, as far as the computation of the degeneracy 
spectrum is concerned, both countings can be formulated in very similar terms as done 
in [|T4|, where they were presented in such a way that the only differences between them 
could be traced back to the form of the projection constraint. 

In section |lll] we have seen that the horizon states in the ABCK framework are de- 
scribed by U(l )-CS labels that can be identified (by using the quantum boundary condi- 
tion and introducing an area interval) with lists of spin components (raplabels) that can 
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be conveniently used to determine the entropy ||2T|. We want to insist at this point on the 
equivalence of the original ABCK entropy definition and the prescription provided in 
| |2T| to compute the entropy by employing only mplabels. Punctures carrying raplabels 
are also characterized by the SU(2) irreducible representation ji associated with the edges 
piercing the horizon. These play a role in the ABCK framework because the compatible 
mi-lists that are counted in the DL prescription must be such that they correspond to lists 
of ji-labels leading to an area eigenvalue in the prescribed area interval. The GM count- 
ing, however, takes the j's as "independent" labels. Although the GM approach has not 
been derived from first principles, according to some authors "[it] is not an unreasonable 
proposal" 111; in fact, we will show that it leads to predictions that are quite close to the 
ones derived by Engle, Noui and Perez in [24] (and in qualitative agreement with the 
ABCK approach). 

If we follow this approach, the prescription to obtain the entropy can be stated as 
follows: 

The entropy S%^{a) of a quantum horizon of the classical area a, according to the GM pre- 
scription, is defined by computing 

S^^(a)=log(1+N^^(a)), 

where N^'^fa) is the number of all the finite, arbitrarily long, sequences ((ji, rai), . . . , (Jn) ttln)) 
of ordered pairs of non-zero half integers ji and spin components mi G {— ji, — ji + 1 , . . . , ji}, 
satisfying 

N 

2 + 1 ) < 0- flwrf ^ mi = 0. 
1=1 I 

Here we want to comment on the type of area interval used in the previous definition. 
Whereas the mathematical details of the DL prescription require that the area interval 
must have the form [0, a], there is no restriction in the GM approach on the form of the 
interval. By invoking the thermodynamical limit, one would expect that the width of the 
area interval is irrelevant. As far as the applicability of our method is concerned, it is 
equally easy to handle an interval of the type [a — 6, a + 5], as originally proposed by 
Ghosh and Mitra, than an interval such as [0, a], because the associated counting can be 
obtained by subtracting the results for [0, a + 5] and [0, a — 6]. 

The combinatorial problems associated with the computation of S<^(a) can be solved 
by essentially following the same steps as in the previous section. These are now: 

Step 1. Fix a value for the area a and obtain all the possible choices for the half integers ji 
compatible with the area, i.e. satisfying 

N 

1=1 

At this point we do not look yet at the ordering of the labels. Notice that this step is 
precisely equivalent to the first step in the DM counting, if one now writes ji = ki/2 
instead of |mi| — ki/2 for the corresponding labels of each configuration c G C(a). 
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Step 2. Count the different ways that we can reorder each of the previous multisets. 

Step 3. For each of the configurations corresponding to the first two steps (ordered se- 
quences of non-zero half-integers), we have to find the possible choices for the mr 
labels that constitute the second element of the pairs that define the configurations 
is this case. Notice that one can think of this step as a generalization of the third step 
in the DL setting. The difference is that there we only had to worry about the signs 
whereas here each mj is a spin component taking values in {— ji, — ji + 1 , . . . , ji}. 

Step 4. Repeat the same procedure for all the area eigenvalues in the interval [0, a] and add 
the number of sequences ((ji, rai)) obtained for each value of the area. 



Step 1: The Pell equation again. 

The first step now is identical to the one described for the DL counting. The only thing 
that has to be remembered is the different meaning of the k/2 labels obtained in the GM 
case. 

Step 2: The reordering degeneracy (r-degeneracy) 

This is again the same as before. The number of different reorderings is the one ob- 
tained in the previous section. An interesting comment that must be made here is that, 
as we will show later, the origin of the non-trivial structure in the black hole entropy and 
the degeneracy spectrum comes from this reordering degeneracy. Hence, we expect to 
have the same type of qualitative behavior both for the DL and the GM counting. 

Step 3: Solving the projection constraint (m-degeneracy) 

Up to this point, we have just determined the possible values of ji = ki/2 in the se- 
quences ((ji,Tai)). In this counting procedure "m-degeneracy" will refer to the number 
of different ways to assign one of the 2ji + 1 values of rai e {— ji, — ji + 1 , . . . , ji} to each 
of the ji labels. This assignment must be subject to the restriction given by the projection 
constraint nij = 0. This problem can be solved in several ways. Some of them are 
analogous to the ones used in the study of the DL counting. However there is an addi- 
tional method -relying on the use of fusion matrices- that suggests an intriguing con- 
nection with conformal field theory and the SU(2) proposals that we will discuss later in 
subsection HV.2I 

The partition problem: The first approach to solving this problem is based on the one 



used in [41] for the resolution of the partition problem as described in III.3 Given the 
sequence (ji , . . . , jN)/ ji G N/2, we want to find the number of sequences (rai) satisfying 
^lii TTLi = M and rai G {—]h ~]i + 1 , • • • > ji}- This number can be computed by intro- 
ducing, as before, an appropriate Hamiltonian for an auxiliary system of spins. In this 

case we just have to substitute the multiples of the Pauli matrices kicrj^' by the operator 
Sf^ = 1(g)l . . -(gjl^sj^'®!®- ■ -(g)!, where si"*' : C^'i+^ C^^i+^ are third-spin component ma- 
trices associated with a jpspin. We consider then the Hamiltonian H : C®'^'^^^' —) C®'^'^^^' 



H = M - ^ Sj 



1=1 



24 



Proceeding as above, we see that 

1 1 T. exp ,-4^m/L, ^ 1 (n(2i. . n . L e-"™- fl ""^^Hl^) 

£=0 \ 1=1 «=1 1=1 / 

{mi=-ji,...,ji} I 

where L is an integer that can be fixed as L = 1 + 2M + ki = 1 + 2M + 2 ji when 
2M + 2 ji is even and L = 2 + 2M. + 2 ji when 2M + 2 ji is odd (in order to 

prevent the vanishing of the denominators). We conclude that the number that we are 
looking for is 

.1=1 f=1 1=1 ^ ' ' 



Fusion matrices: The second approach relies on techniques used in the context of con- 
formal field theories | |42| and, in fact, is suggestive of a deep connection between them 
and the problem of computing black hole entropy in LQG ||25|. The key insight now 
is to realize that the problem of determining the number of solutions to the projection 
constraint ^jTTLi = is equivalent to counting the number of SU(2) irreducible repre- 
sentations, including multiplicities, that appear in the tensor product In the 
following we will denote the SU(2) irreducible representation corresponding to spin ji as 
[ji] = [ki/2]. 

We start by writing the tensor product of two SU(2) representations as 



ki 
2 



2 



k3=0 



2 



where the integers ^NfJ^^,^^, called fusion numbers [|42|, give us the number of times that 
the representation [k3/2] appears in the tensor product of [ki/2] and [k2/2]. For each 
k G N U {0}, we introduce then the infinite /uszon matrices (Ck)k,k2 
k2 G N U {0}. These satisfy the recursion relation 

Ck, k = 0,l,... 



:= KJ^^^, where ki. 



C 



k+2 



xc 



k+1 



(IV.l) 

where we use the notation X := C^. Explicitly, Xk,^^ = 6ki,k2-i + Sk,,k2+i/ which shows that 
X is a so-called infinite Toeplitz matrix II43B . The solution to (IV.l I, with initial conditions 
Co = I and Ci = X, is 

Ck = UjX/2), k = 0,l,... 

where the are the Chebyshev polynomials of the second kind. The tensor product of 
any number of representations can be decomposed as a direct sum of irreducible repre- 
sentations by multiplying the fusion matrices defined above. This way we get 



ki" 




'k2" 




'kN" 








8) ■ ■ ■ (g) 


_2_ 




_2_ 


_ 2 _ 



k^ ■ ■ ■ )ki k 



k=0 
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Notice that the product of matrices 



Ck2 



Ck, = Uk, (X/2)Uk3 (X/2) • • • Uk, (X/2) 



is a polynomial in X. The total number of representations giving the solution to the com- 
binatorial problem that we are trying to solve now is 



y~ i C]c2 Ci< 



CkN)k,k 



2^(Uk,(X/2)Uk3(X/2) ■ • • UkjX/2))k,k, (IV.2) 



k=0 



k=0 



i.e. the sum of the (finite number of non zero) elements in the ki row of the infinite matrix 
Ck2 Ckj ■ ■ ■ CkN . An integral representation for this sum can be obtained by introducing, as 
in [|43|, a resolution of the identity for the Toeplitz matrix X and the identity Uk(cos 9) = 
sin[(k + 1)9]/ sin 9 for the Chebyshev polynomials. Equation (IV.2) can be equivalently 
written as 



k=0 



9 



9 



^(Ck^Ckj ■ ■ ■ CkN)k,k = - cos - fcos--cos(^ki + -)9 j Y[ 



1=1 



1=1 



sin(ki + 1)9 
sin 9 



d9. 

(IV.3) 



Generating functions: The third way to solve the projection constraint for the GM 
counting makes use of generating functions. The idea is to take now 



n L 



,k-2n 



n 



.n=0 



,k+l 



r-k-1 



This function is similar to the one corresponding to the DL counting ( |III.16 >. The number 
of solutions to the projection constraint taj = M is given now by the coefficient of z^'^ 
in the previous expression and can be extracted as before by using Cauchy's theorem. 
The relevant coefficients for a given 2M are given now by 



[z^"^] n 



,k+l 



,-k-l 



1 

27ti 



() 



dz 



n 



(IV.4) 



where C is, again, an index-one curve surrounding the origin. By choosing for it a unit 
circumference, the previous integral can be written in the following useful alternative 
form 



[z'^] n 



,k+l 



,-k-l \ Nk 



r-1 



2n 



-TT sin^''(k+ 1)9 ,^ 
cos(2M9) _Q \, ^ d9 



sin'^^ 9 



(IV.5) 



In Appendix |C] we give yet another procedure to solve the projection constraint based on 
group theoretical methods [see equation (C.l) and compare equations (IV.5 1 and (C.2|]. 



GM black hole degeneracy spectrum. 
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As we did in the DL case, given an area a = ^ • qty^ we can define the m-degeneracy 
d^;;'^(c) of the configuratiori|^c G e(a) as 



sin^^i^ (27t(ki,+ 1)£/L) 



e=i i a 



sin ^ix(27r£/L) 



2 rn 
71 





1 



cos^ cos;-cos(^^kiN,,+f)e nn 



sin^^Kk^„ + 1)9 
Sin ^cc 9 



de 



In 



prpr sin^^Kk-,+ 1)9 
i.^^i sin^^^.9 



where the integer L = L(c) has been defined above. 

We will also define the black hole degeneracy associated with a given value of the area 
a = Li qiVP^as 



cee(a) 



where d.r(c) was given in (III.9). In the following we will refer to D*^ (a) as the black 
hole degeneracy spectrum for the GM counting. According to the previous discussion, for 
every value of the area a, it gives the number of sequences ((ji, rai)), where ji G N/2 and 
mi G Z/2, satisfying the conditions 



^ Y VUiU + 1) = ci) TUi G {ji -n : n = 0, ...,2ji}, ^ tui = . 



1=1 



1=1 



Step 4: Adding up. 

In the previous paragraphs we have given a procedure to compute the black hole de- 
generacy spectrum for the GM counting. As in the DL counting, we have to add the 
degeneracies corresponding to all the possible values of the area spectrum in a certain 
interval and, again, use integral transform techniques to deal with the difficult combina- 
torial problem of summing for all the relevant values of the area. To this end we proceed 
as in subsection |III.4[ Let us consider the generating function for the black hole degener- 
acy D^'^(a) (see Appendix|A]) 



-1 



G^^(z,xi,x„...)= jx^M . (IV.6) 

\ 1=1 m=l / 

and perform the substitutions z = e^^, xi = e^^^, xi = e^^^ , . . . , to get 

'l - f ^M^+^)^ ,-sVm^] '\as d3 do. . (IV.7) 



expS^^(a) ^ 



rln 


fXo+ioo / 






. . 


Xq— too \ 



sm cu 

k=l 



^ The definition of the set of configurations 6 (a) here is the same as in the DL scheme and relies on the 
solutions to the Pell equations as before. 
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This expression corresponds to an area interval of the form [0, a] in the definition of the 
entropy. As explained before, the values for an area interval [a — 5, a + 5] can be easily 
obtained from it. The behavior of S^'^fa) as a function of area is the same as for the DL 
counting, namely a linear dependence and a logarithmic correction with —1/2 coefficient. 

The Bekenstein-Hawking entropy-area relation is recovered now by taking y = 0.274 

As in the case of the DL counting, it is possible to define a simplified combinatorial 
problem by ignoring the projection constraint and considering all the 2j + 1 choices of 
m-labels |[T9ll22| . This can be easily done just by replacing sin(k + 1 )a)/ sin cu by k + 1 in 
equation ( |IV.7| ). 



IV.2. SU{2] black hole entropy 

The initial proposals of Smolin and Krasnov ||28l |30| to study black hole entropy in 
LQG (predating the more detailed ABCK approach) suggest that the horizon degrees of 
freedom could be accounted for by a SU(2)-CS theory. One of the non-trivial points of the 
ABCK treatment is the claim that the horizon degrees of freedom are, in fact, described 
by a U(1 )-CS model. This has been justified in this framework by a combination of ideas 
involving the analysis of gauge transformations on the isolated horizon and Hamiltonian 
methods. However, the definition of an isolated horizon does not require the introduction 
of any internal symmetry so, on the face of it, no restriction on the internal symmetry on 
the horizon is expected when connection variables are used [24] . The subtleties associated 
with this issue have led some authors to explore the possible consequences of ascribing 
the black hole degrees of freedom (and, hence, the entropy) to a SU(2)-CS theory. In 
particular, Kaul and Majumdar [|23l |33| have argued that the Bekenstein-Hawking area 
law can be explained in a such scheme and have also proposed logarithmic corrections 
to the black hole entropy. The entropy computations carried out in [23, 33] rely on well- 
established conformal field theory methods [42] and, hence, will not be reviewed here. 
The main problem with this approach is that the SU(2)-CS model on the horizon is really 
put in by hand and not derived from the quantization of the relevant sector of general 
relativity. Although it is definitely interesting to explore alternative approaches it stands 
to reason that any model aspiring to provide a fundamental physical description of black 
holes should be solidly rooted in theory. 

Other authors have attempted to see if a SU(2)-CS description can be derived from 
first principles within LQG. In particular Engle, Noui and Perez (ENP) have proposed a 
covariant Hamiltonian scheme leading to a pure SU(2) formulation [2^ (see also [34]). 
These results are not compatible with the ABCK model but lead to a remarkable entropy 
definition in the sense that it gives the same Immirzi parameter predicted by Ghosh and 
Mitra and the same logarithmic correction proposed by Kaul and Majumdar. 

The entropy in the ENP model is obtained according to the following prescription: 

The entropy S™^(a) of a quantum horizon of the classical area = Arcyl^K (when y < \/3) 
is defined as 

Sr(a)=log(1+NfP(a)), 
where N™''(a) is the number of all the finite, arbitrarily long, sequences [ji, . . . , jN) of non-zero 
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half integers ji satisfying 

N 

1=1 

and counted with a multiplicity given by the dimension of the invariant subspace Inv(®i[ji]). 

The last condition in the ENP entropy prescription plays the role of the projection con- 
straint characteristic of the other approaches. Actually, if instead of using Inv ( ® i [j i] ) each 
sequence is counted with a multiplicity given by the number of irreducible represen- 
tations -taking into account multiplicities- that appear in the direct sum decomposition 
of the tensor product ®i[ji] we recover precisely the GM counting. Hence, the number of 
configurations accounting for the entropy for a given area in the ENP approach is slightly 
lower than the one in the GM counting. It only changes the logarithmic correction (see 
below). 

We discuss now in an schematic way how the ENP entropy is computed in our scheme. 

Step 1. Fix a value for the area a and obtain all the possible choices (without ordering) for 
the half integers ji compatible with the area, i.e. satisfying 

N 

^Viiih + ^) = j. 

1=1 



Step 2. Count the different ways to reorder each of the previous multisets. 

Step 3. Determine the dimension of the invariant subspace Inv(®i[ji]) associated with each 
sequence (ji). 

Step 4. Repeat the same procedure for all the area eigenvalues smaller than a and add the 
numbers obtained in each case. 

The only difference with the GM case is in step 3. In prac tice f4^, this means that it 
suffices to replace the m-degeneracy d^'^(c) by (see equation (C.3| in Appendix|c]) 



1 



27ti 

1 r^TT 



dz (z — z 
c ^ 



-1 ^2 



() 



2 



nn 
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sm 







We will also define the black hole degeneracy associated with a given value of the area 
^ = Liqi^Pi as 

cee(a) 

where dr{c] is given by (III.9). We will refer to D™^( a) as the black hole degeneracy spectrum 
as before. The entropy can be computed by summing up black hole degeneracies. In 
practice, this can be done by considering the generating function [|44j 



'ENP, 



Z,Xi,X2, 



2 



i=l Ta=l 



r-ki„-l 



-1 



.-1 



(IV.8) 
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and performing the substitutions z = e^^, xi = e xi = e . . . , to get 

c2n rXfl+ioo 



exp S 



ENP/ 



a 



[InYi J 



s sin^ UJ 



- )- ^^^^+^^^ e-v^)'^e-dsdcu. 
V — sincu / 



xq-ioo i^^I 

This expression leads us to fix y = 0.274 . . . and gives the logarithmic correction 

3 



2 



log(a/£^p). 



In Figs. |3|and|4|we plot D^^^fa) and S<^^ ( a) as functions of the area. Notice that these 
figures are similar to Figs. [T]and|2] 



D™''[a} 




FIG. 3. Plot of the black hole degeneracy spectrum D (in units of 1 ) in terms of the area (in 
units of 47ty£p) for a range of area values. 



V. SUMMARY OF THE DIFFERENT SCHEMES 

The behavior of the entropy as a function of the area, obtained by using the differ- 
ent countings considered in the previous sections, is compatible with the Bekenstein- 
Hawking law for appropriate choices of the Immirzi parameter. The qualitative behavior 
for the different models is shown in Figs. |2] and |4] (notice that the GM and ENP countings 
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a 



FIG. 4. Plot of S™^(a) in terms of the area (in units of 47ty£p). This plot is essentially the same for 



differ only in the logarithmic corrections and, hence, the plots of the entropy in these two 
cases are essentially indistinguishable). These figures show that the entropy grows in 
steps of a characteristic width Aa that appears to be the same for all the countings (when 
measured in units of 47ty£p), i.e. 



2^ 
4n 



(Aa 



DL 



(Aa) 



GM 



[Aa] 



ENP 



8.78- 



47ryDL«^ 47ryGM«p 



47ryENp£p 



4n 



0.69 ■■ 



There are also subdominant logarithmic corrections that differ from model to model. 
These results are summarized in the following table. 



Approach 


y 


Logarithmic correction 


Aa 


ABCK-DL 


yoL = 0.237- ■ ■ 


-1/21og(a/£i) 


XYdl^p 


GM 


ycM = 0.274- ■■ 


-1/21og(a/£2) 


XTgm^p 


ENP 


Tenp = Tgm 


-3/21og(a/£2) 


XTenp^p 



A comment regarding Aa is in order now. In the previous sections, we have measured 
areas in units of Anylp. This choice of units is, in principle, model-dependent because 



31 



the value of y must be chosen differently in order to arrive at the Bekenstein-Hawking 
law. This means that area steps in the entropy have different (but related) sizes for the 
different countings and hence the actual value of y leaves a characteristic signature in the 
behavior of the entropy A final remark regarding x is the conjecture appearing in [,20l 
stating that x = 8 log 3. 



VI. DETAILED ANALYSIS OF THE ENTROPY MICROSTRUCTURE 

The main goal of this section is to use the techniques presented in the preceding parts 
of the paper to obtain a detailed understanding of the behavior of the entropy The "mod- 
ular" nature of the number- theoretic approach that we have followed will allow us to 
ascribe some of the interesting behaviors that we want to understand (in particular the 
staircase growth of the entropy) to some of the specific steps that we have followed. This 
analysis is valuable on its own but can provide, in addition, important information to 
tackle the difficult problem of fully understanding the asymptotic behavior of the en- 
tropy in the large area limit. 

We are going to separate the analysis in different parts. On the one hand, we are 
going to study the relative importance of each of the possible values of the k-labels in 
the different configurations. On the other hand, we are going to analyze the contribution 
of each type of degeneracy, namely the r-degeneracy and the m-degeneracy, and their 
influence in the observed patterns. 

From a practical point of view it is better to focus on the behavior of the degeneracy 
spectrum D(a) instead of the entropy. In the following we will concentrate on the DL- 
counting because the arguments can be copied for all the other models that we have 
discussed in the paper. 



VI.l. Disentangling the area spectrum with the help of squarefree integers. 



The methods developed in section III will allow us to analyze in detail the structure 
of the black hole degeneracy spectrum and the entropy. The area spectrum is given by 
certain integer linear combinations of square roots of squarefree numbers a = '^iy/Vi- 



The solutions to the diophantine equations described above (subsection III.1[ > tell us how 



to identify the configurations C(a) associated with each value of the area. One of the 
important features of our formalism is the central role of the squarefree numbers in the 
classification of the points in the area spectrum. In many instances, it is much more 
convenient to use them instead of spin labels. Of course, once a squarefree is fixed, there 
is an infinite number of possible spin labels given by the solutions to the Pell equation 
(III.7I. The squarefree numbers and the associated spin labels k/2 are given in table|lj The 



i label appearing in the first column identifies the elements of the sequence of squarefree 
numbers listed in increasing order (pi = 2, p2 = 3, ps = 5, . . .). The squarefree numbers 
themselves are shown in the second column. They are listed according to their order 
of appearance in the successive eigenvalues of the area spectrum. The spin labels in 
the column are associated with the fundamental solution of the corresponding Pell 
equation. The spins in this column grow in units of 1 /2. This pattern is only interrupted 
by the appearance of "secondary" solutions in k2, kj, . . . (for example k/2 = 3 appears as 
the second solution to the Pell equation for p2 = 3, i.e. in the ^ column). 
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As we know form previous analyses [ jSTl |22| . the smaller the k-label is, the higher its 
contribution to the most degenerate states. This is the reason why for many practical pur- 
poses we can concentrate on the study of the lowest squarefree numbers. For example, 
if we want to consider only spin labels up to 5/2, we only have to take the squarefree 
numbers 3, 2, 1 5, 6, and 35. The fact that the very first value of k, (corresponding to V3) is 
1 /2 justifies the well known observation that punctures with spin labels equal to 1 /2 give 
the most significant contribution to the entropy. 

Subfamilies in the area spectrum consisting of area eigenvalues that are multiples of 
the square root of a single square free number are specially easy to analyze. Figure 
[s] shows the behavior of logDj!^^(a) and \ogD^^[a) in three different cases given by the 
subfamilies of the area spectrum that can be written as integer multiples of V5, V2 and 
VTS. As we can see each of these subfamilies is characterized by a different growth rate. 
The logarithm of the black hole degeneracy spectrum grows linearly with the area with 
a slope Si that can be determined numerically for each squarefree number by finding 
the unique real solution to the equation |[T8| 

oo 

1 _2^e"''v^^-"'-+^' = (VI.l) 

a=l 

in the unknown s^. The corresponding values are given in the last column of table |l| As 
we can see they decrease monotonically for each successive pt. In addition to the linear 
growth there is a logarithmic correction of the form — (log a)/2 in the physically relevant 



TABLE I. Pell equation: area spectrum and squarefrees 



t 




\c\/2 


ki/2 


kl/2 




Si 


2 




1/2 


3 


25/2 




0.4526794 


1 


V2 


1 


8 


49 




0.2545910 


10 


\/T5 


3/2 


15 


243/2 




0.1808782 


4 


V6 


2 


24 


242 




0.1418796 


22 


V35 


5/2 


35 


845/2 




0.1172453 


5 


V7 


7/2 


63 


2023/2 




0.0873322 


3 


V5 


4 


80 


1444 




0.0774971 


7 


\/TT 


9/2 


99 


3969/2 




0.0696641 


18 


Vso 


5 


120 


2645 




0.0632754 


88 




11/2 


143 


6875/2 




0.0579639 


27 




6 


168 


4374 




0.0534775 


119 


V^95 


13/2 


195 


10933/2 




0.0496373 


9 




7 


224 


6727 




0.0463128 


156 


V255 


15/2 


255 


16335/2 




0.0434066 


198 


V323 


17/2 


323 


23273/2 




0.0385677 
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case for which the projection constraint is taken into account. It is important to point out 
here that the actual growth rate of each family can only be determined once the value of 
the Immirzi parameter y has been fixed by enforcing the Bekenstein-Hawking area law. 

The asymptotic approximation for the subfamilies considered above can be easily ob- 
tained from the generating functions discussed in Appendix A We will describe next how 
this can be done for more general subfamilies of points in the area spectrum consisting 
of integer linear combinations of a fixed number of square roots of squarefree numbers 
Vu Pin labeled by the subset J = {ii , . . . , in} C N. The idea is to restrict the full generat- 
ing function ( IIL24[ > by eliminating all the variables with the exception of a finite number 



of them, (xi, , . . . ,x^^), associated with the square free numbers defining the subfamily. 
This way we get 

oo ^ 



ieJ <x=^ 



By using the Laplace transform methods described in flSl we can easily find the asymp- 
totics corresponding to the subfamilies by performing the substitutions xi = e^^^, z — 1 , 
and solving the equation 



^ -2^Y_ = (VL2) 

tea a=i 

that gives the real pole Sj. This directly provides the slope of the asymptotic linear growth 
of the degeneracy spectrum associated with the chosen subfamily. The functions of the 
type 

oo 

iea a=i 

have some distinctive properties that justify the previous procedure: 

• They have a single real zero sj. 

• The remaining complex zeros have real parts that are smaller or equal than sj. 

• The real parts of these zeros accumulate to the value Sj. 



These properties explain why the real solution to equation (VL2) captures the effective 
asymptotic behavior of the black hole degeneracy for the chosen subfamily, and leaves 
room for interesting (i.e. periodic) substructures. Several comments are in order now: 

• The derived growth applies to the "convex" envelope associated with the chosen 
subfamily. 

• There may be additional substructures (in the form of bands for the degeneracy 
spectrum or steps for the full entropy). 

• There may be subdominant (i.e. logarithmic) terms in the asymptotic expansions. 
To find these the z variable must be reintroduced. 



For a given subfamily there are sub-subfamilies (chosen by further restricting the 
number of squarefrees) with different growth rates. 
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As we have mentioned above, the behavior of the black hole degeneracy spectrum for 
area eigenvalues that can be written as integer multiples of a single squarefree is shown 
in Fig. |5] The expected linear growth (with the corresponding logarithmic corrections 
when the projection constraint is included) can be readily seen. 

Figure [6] shows the degeneracy spectrum for the subset of area eigenvalues consisting 
of integer linear combinations of ^Jl and V3 when the projection constraint is included. 
There are several interesting features that can be seen in this case: 

• The set of points represented in Fig. [6] corresponds to a subset of the area eigenval- 
ues consisting in the union of subfamilies ^i(q°) of the form 

^i(q2)={qiV2+q2V^: q2GNU{0}}, A^[c^\^ = {^\^ ^ c^2^ : qz G N U {0}} 
where q? are fixed non-negative integers. 

• The asymptotic growth of the degeneracy spectrum within each family .Ai(q°) is 
controlled by Si. 

• Notice, however, that the "envelope" of the plot grows with a rate given by the 
Sj = 0.645008... with J = {1,2}. This is larger than both si and S2. Also notice 
that some points in each of these families contribute to the envelope that eventually 
defines the growth of the degeneracy spectrum for intermediate values of the area. 

• Although the band structure of the full degeneracy spectrum is not apparent at this 
stage, the separation in area of the points that define the envelope is quite close to 
the actual periodicity of the full spectrum. In fact it is possible to argue that these 
points somehow define the "roots" of the bands. 

Finally, Fig. |7] shows the result of considering area eigenvalues that can be written as 
integer linear combinations of \fl., and vT^. The appearance of the band structure in 
the degeneracy spectrum is evident now. Figure |8] shows a detail of the plot for larger 
values of the area spectrum where the band structure can be clearly seen. The envelope 
of the graph grows now with a slope given by Sj = 87 ■ ■ ■ which is actually quite 

close to the actual value obtained for the full area spectrum sn = 0.746232 ■ ■ ■ Several 
comments are in order now: 

The first is to notice that the band structure in the spectrum appears clearly for area 
eigenvalues involving a small number of square roots of squarefree integers. In fact, with 
linear combinations of the first two of them {^/3 and \/2) it is not possible to distinguish 
those structures in the spectrum, but as soon as the first three are taken into account, the 
bands readily show up as can be seen in Figs. [7|and|8j 

The second point to make is that, once the first three squarefrees are considered, 
adding more linear combinations does not significantly change the structure and position 
of the bands. As one adds new squarefree numbers to the values of the area spectrum, the 
only effect that can be observed is that more and more points appear 'in' the peaks, but 
not 'out' the peaks, giving as a result a better resolution of the shape, but not modifying 
any feature of the structure in a significant way. 

The third point is that, as soon as the five or six first squarefree numbers (in the order 
given in table |l]) are considered, the spectrum obtained is almost indistinguishable from 
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the complete one. This means that, with the first five squarefree numbers, we are able 
to reproduce all the points in the spectrum that have a degeneracy of the same order of 
magnitude as the maximum. A possible way to quantify this fact is by looking at the 
asymptotic growth rate of the envelope given by ( |VI.2 > as a function of the number of 
squarefrees used. This is shown in table |Il| The information that one can extract from 



TABLE II. Asymptotic growth rate 



S2 


0.452679 . . . 




0.645008 . . . 


S{1,2,10} 


0.708187... 


S{1,2,4,10} 


0.732215... 


S{1,2,4,10,22} 


0.742270 . . . 


5(1,2,4,5,10,22} 


0.744388 . . . 


5(1,2,3,4,5,10,22} 


0.745368 . . . 






Sn 


0.746232 . . . 



the above analysis is valuable, in particular the fact that the highest degeneracy states 
responsible for the observed pattern in the degeneracy spectrum are (mostly) composed 
by punctures with k = 1,2,3,4,5, and 6. Moreover, we know that the area values cor- 
responding to these configurations are those containing multiples of square roots of the 
squarefree numbers 3,2, 15, 6, and 35. Then, considering only these values of the area, 
we are sure that no relevant information is being missed but the calculations are simpli- 
fied. In fact, the plots appearing in the present paper corresponding to the largest area 
eigenvalues considered have been obtained by using this approximation. 



VI.2. Effects of the r-degeneracy and ra-degeneracy on the entropy microstructure 

Irrespective of the model (DL, GM, ENP), the black hole degeneracy spectrum D ( a) 
has the form 

D(a) = Y_ cir(c)d^(c). 

cee(Q] 

There are then two different contributions to the degeneracy of each configuration: one, 
the r-degeneracy, is model independent and comes from the possible reorderings of its 
labels. The other, the ra-degeneracy, originates in the projection constraint (or similar 
conditions). In this subsection we want to study the relative contributions of each of 
them and, specifically, try to understand if any of these can explain by itself the observed 
microstructure of D(a). To this end we will look at some auxiliary objects built out of the 
dr(c) and drn,(c), in particular 

D,(a) := d,(c) , D^(a) := dn,(c) . 

cee(Q) cee(Q) 

Figure |9] shows Drafa) for the DL case (the other cases behave in a similar way). As 
we can see, there is some non-trivial substructure in this plot but this does not resemble 
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FIG. 5. Plot of log DP^(a) and log D^^(a) in terms of the area (in units of Anyip). The figure shows 
the values associated with area eigenvalues multiples of Vi, \f2. and \/T5 respectively. The solid 
lines correspond to the asymptotic approximations given by straight lines with the slopes S2, si 
and sio appearing in table |l] When the projection constraint is taken into account (lower plot) the 
asymptotic approximations have also a logarithmic correction — (log a)/2. 
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FIG. 6. Plot of log D^^(a) for area eigenvalues of the form qi \/2 + q2\/2 with qi , qi G N U {0}. The 
solid line with the largest slope S{i 2} = 0.645008 • • • (determined by equation ( VI.2 1 for J = {1 , 2}) 



gives the asymptotic approximation for their growth. The other solid lines correspond to the 
\/2 and \/3 subfamilies. As it can be seen the "envelope" of the points grows faster than these 
subfamilies. In all the cases the logarithmic correction — (log a) /2 has been included. 
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FIG. 7. Plot of the logarithm of the black hole degeneracy spectrum, logD^^(a), in terms of the 
area (in units of Anylp). The figure shows the values associated with area eigenvalues that can be 
written as integer linear combinations of VS, \/2 and VTS. The arrows mark the position of the 
bands as predicted by equation ( VI.4| , taking into account that only even values of P contribute 
when the projection constraint is considered. Notice that beyond a = 5 there are bands at each of 
these positions. Notice also that there are many values of the area spectrum for which D^^(a) = 
and, hence, do not contribute to the band structure shown here. The colors indicate the subfami- 
lies to which the different points in the plot belong (blue for \/3, red for \/2, green for \/T5/ purple 
for linear combinations of ^/2. and VS and black for linear combinations of \/2, \/3 and \/T5). 



the band structure that can be seen in the full degeneracy spectrum D(a). On the other 
hand, the plot of Dria] (Fig 10 1 readily show^s a band structure similar to the one found 
in D(a). The only significant difference is that the exponential growth of the peaks is less 
pronounced. Therefore, we can conclude that the origin of the band structure resides in 
the r-degeneracy, which is rooted in the distinguishable character of the punctures. This 
is a consequence of the action of general diffeomorphisms over the horizon states. The 
fact that this degeneracy is common for all the countings explains why the band structure 
appears for all of them. 



VI.3. The peaks in the degeneracy spectrum 

In order to study the appearance of the peaks in the degeneracy spectrum it is very 
useful to introduce the following function in the space of black hole configurations ||20| 



P(c) :=3K[c) + 2N(c), 
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FIG. 8. Plot of (in units of 1 0^°) versus area (in units of Anylp) considering only areas involv- 
ing the squarefrees numbers 2, 3, and 15. 



where, given c = {(Tc, N]^)}, we define K(c) := ^^kN^ and N(c) := ^^N]^ (the sums 
extend to all the k-labels appearing in c). P(c) is always a natural number, however, it is 
important to notice that those configurations giving rise to odd values of P(c) correspond 
to zero degeneracy. This is so because the projection constraint cannot be satisfied if P is 
odd. 

As we justify now, the crucial property of P(c) is that each peak in the degeneracy 
spectrum is characterized by a single value of P (see Fig. [Tl]). This fact can be under- 
stood by following a simplified analysis, that we explain next, that captures the essential 
features of the problem. To begin with let us consider the function defined by 



Ff(c) :=d,(c)d°Hc) = 



in the configuration space C = C(a). Here we are neglecting the projection constraint 
(though we will restrict ourselves t o cons idering only even values of P). Notice, by the 
way, that the degeneracy given by (IIL20I is just the sum of FjJ^^fc) for all the configura- 
tions c G C(a) C C corresponding to a specific area a. Now, it is possible to partition 
the space of configurations according to the values of P(c). This procedure should be 
thought of as a coarse grained alternative to the standard approach of describing the de- 
generacy spectrum in terms of a because, for a fixed value of P, there are configurations 
corresponding to a range of areas. The remarkable fact, as we will see, is that -in the 
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FIG. 9. Plot of B^{a) := Lcee{a) dm^(c) (in units of lO") versus area (in units of 47Ty£^) 



continuous approximation- there is a single maximum of F^'^fc) for each value of P. Fur- 
thermore, the areas associated with the configurations that maximize F°^ for the different 
values of P closely match the observed position of the peaks in the degeneracy spectrum. 

We show why this is so by considering continuous values for the variables Ni^ and 
extending the function accordingly. We then use the Lagrange multipliers method to 
enforce the condition that P(c) = P with P G N. By taking logF^'^fc), using the Stirling 
approximation, introducing a single Lagrange multiplier A, and using the notation tslk := 
Ni^/(^^, Nv), the extrema are determined by 

ISlk = 2exp(-A(3k + 2)) . (VI.3) 

Adding the previous expressions in the label k we get the condition 

l=^2exp(-A(3k + 2))=^^ 



that fixes the value of A = 
quintic equation 



— log yo in terms of the single real solution of an auxiliary 
2^5 + - 1 = . 
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FIG. 10. Plotof r-degeneracy Dr(a) = Hcee(Q) d^rtcjSTatc) (in units of 10^^) versus area a expressed 
in units of 47ty£p. Here 6m (c) is one if the configuration c satisfies the projection constraint and 
zero otherwise. Without this factor the period of the peaks decreases by a factor of two. 



The solution Vq can be written in terms of hypergeometric functions as 

2 , 



4F3 
16 



1 4 7 13_ 1 2 7_ 6250\ 
15' 15' T5' 15' 3' 3' 6' iT J 



4r3 



2 3 4 6 2 4 3 _6250\ 
5' 5' 5' 5' 3' 3' 2"' 2^7 



+ y4F3 



11 14 17 23 4 5 11 _6250\ 
T5' 15' 15' 15' 3' 3' T' 27~ ) 



and it gives A = 0.260847 ■ ■ ■ 

Introducing back the value of A in (VI.3I, we obtain the distribution for the fvlk- This 
determines the values of the -maximizing log F^'^ for a given P- modulo the value of 
the sum Nk that can be obtained by using the constraint 

00 00 00 y-, 

P = 3K(c) + 2N (c) = ( }^ Nk) ■ (2 + 6 }^ ke-^(3'^+^') ^ }^ = ^ 00 

k=l k=l k=l 2-k=] - 

Finally, the area corresponding to the maximum of in terms of the peak label P is 



-A(3k+2) 



a(P) 



k=l 



,-A(3k+2) 



1 + 3 1 ke-^(3k+2) 



P = (0.34959- -OP. 



(VI.4) 
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Several comments are in order now. First, it is interesting to compare the distribution 



given by (VL3) with the so called "maximum degeneracy distribution" (MDD) given by 
f^k = 2exp ( - ADL\/lc(k + 2)) , where Adl = 0.746232 . . . (VL5) 
The MDD can be obtained by maximizing log F^'^ subject to the condition that the area 



+ is fixed. As can be seen in Fig. 12l they are very close to each other 



and, in particular, they are almost identical for the most relevant spin labels. A second 



remark concerns the positions of the a(P). As can be readily seen in ( VI.4| they are integer 
multiples of the constant 0.34959 ■ ■ • , and in the physically relevant case where the pro- 
jection constraint is taken into account (even P's), they are multiples of 0.699\ 8 ■ ■ ■ These 
positions are very well correlated with the observed maxima of the degeneracy distribu- 
tion. It is interesting to remark here that a very similar distribution has been obtained by 
a completely different method in [20J. This compatibility is a non trivial cross check of 
both approaches. Third, the previous computations have been performed in a continuum 
approximation that neglects the fact that the values of are integers. In practice, this 
means that the actual shape of the peaks may be more irregular and secondary maxima 
may be present. 

Once we have checked that P plays an important role in the identification and labeling 
of the peaks appearing in the degeneracy spectrum it is straightforward to write down 
generating functions that allow us to select the contribution to the degeneracy spectrum 
of the configurations associated with a given value of P. Following the approach de- 
scribed in Appendix A it is possible to take care of the values of P by introducing a new 
variable y in the generating functions. This way we find 

/ oo oo \ 

[v; z,x^,x^,...)={^-Y_Y_ ^^^''^^[^^^ + z-^'^)xf . (VI.6) 



i=l K=l 



The coefficient 



D°L(a,P) := [z°][x?^X2^^---][vnG^^(A.;z,xi,X2, 



gives the sum of the degenerations d.r(c)d5;; (c) of all the configurations with area a = 
Y-i 1i satisfying P(c) — P. Obviously, the black hole degeneracy is given by 

oo 

D°L(a) =_^D°L(a,P). 

p=i 

A graphical confirmation of the ability of the P(c) to identify the peaks in the black hole 
degeneracy spectrum is given in Fig. qT| Here we have represented D°'-(a, P) as a func- 
tion of a for several consecutive values of P. 

The preceding analysis can be easily carried out for the other countings presented in 
the paper with similar results. As we have shown the band structure is explained by the 
r-degeneracy. The m-degeneracy only plays the role of suppressing some configurations. 
The main consequence of this is the effective change of the periodicity observed in the 
black hole degeneracy spectrum. 
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FIG. 11. This is a new version of Fig. [Tjwhere we have higWighted one of the peaks singled out 
by the peak counter P. We also show the position of the peaks according to formula ( VI.4[ |. 



VII. CONCLUSIONS 



We have given a comprehensive account of a number of combinatorial methods that 
allow us to study in great detail the behavior of the entropy in loop quantum gravity. 
The power of the techniques that we have employed can be appreciated by using them 
to study the different counting proposals that appear in the literature. The number the- 
oretic methods based in the solution of diophantine equations provide a very efficient 
way to study the spectrum of the area operator and characterize the configurations of 
black holes that determine the entropy These methods also provide powerful classifying 
criteria to disentangle the black hole degeneracy spectrum that can be exploited to under- 
stand the origin of the entropy microstructure. This detailed information can be codified 
in an extremely efficient way in suitable generating functions from which it is possible 
to extract detailed information at will. We illustrate the power of these generating func- 
tions by providing complete computations for the smallest black holes in Appendix |B| 
These computations show in a convincing manner how the interesting microstructure of 
the entropy comes to life and show that it is a robust feature present in all the different 
schemes. We give also in appendix|C]a new method based in group-theoretic arguments 
that can be used to obtain, in a unified way, the results concerning the implementation of 
the projection constraint. 

The final goal of the approach that we have described here is to verify if the entropy 
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k 



FIG. 12. Plot of the degeneracy distribution given by formula ( VI.3| (dashed line) compared to the 
maximum degeneracy distribution ( VI.5| (solid line). 



discretization holds for macroscopic black holes. An important tool in this analysis is the 
generating function given by VI.6 because it isolates the configurations that contribute to 
the bands that appear in the black hole degeneracy spectrum. Actually the proper iden- 
tification of this generating function is one of the main new results of the present paper. 
We want to mention that, although this function is written here for the DL counting, the 
method described in appendix A can be used to derive similar expressions for the other 
countings mentioned in the paper. The main open problem -that we expect to solve 
with the methods described here- is the persistence of the observed entropy structure for 
macroscopic black holes. If this is the case, then LQG would realize in a very non-trivial 
way the predictions about the effective equally spaced quantization of black hole areas. 
In our opinion this would lend an important support to the theory. 
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Appendix A: Generating functions 



In this appendix we give a pedagogical account of how the different generating func- 
tions mentioned in the paper are derived (a general treatment of generating functions 
can be found in Il45l ). Our construction uses the set of allowed configurations C(a) con- 
sisting of all multisets c = {(kj^, Nkj^)} associated with a value of the area a = (]i^/pi 
(as explained in section III.ll. The black hole degeneracy D(a) can be obtained from 
C(a) by incorporating two sources of degeneracy The first one is related to the possi- 
ble reordertngs (r-degeneracy) in each multiset and is common for all the models. The 
second corresponds to the additional conditions that define the different countings (ra- 
degeneracy). 



A.l. Generating function for the r-degeneracy 

Given a value of the area a = qi^Pi the r-degeneracy is given by 



Dr(a) = Y_ 



cee(a) 



(A.1) 



where Ni<.i are the solutions to (IIL8 1. In order to find the generating function Gr for Dr(a) 



we split the problem in three steps: 

1. Find the generating function for the number of solutions of each of the diophantine 
equations (III.8), that is, the number 



cee(a) 



2. Modify this generating function to introduce the denominators appearing in the 
definition (|A.1| of D^(a), that is, the number 



L 

cee(Q] 



1 



3. Modify the generating function obtained in the previous step to account also for the 
numerators in ( A.l| >. 



The first step is solved by getting a generating function that counts the number of 
solutions to diophantine equations of the form 



(A.2) 



where the unknowns are the Ni^ and the y\ are given by the solutions to the Pell equa- 
tions. If we formally consider the function 



F^(x) = (x°-yi + x^-yi + x^-y' + ■■■)■ (x°-y^ + x^-y^ + x^-y^ + ■ ■ ■ ) 
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it is not difficult to see that the coefficient of the x'^' term of the McLaurin expansion of 
Fi(x) gives precisely the number of solutions of the diophantine equation ( A.2 1. Now it is 
possible to write F^(x) as 



h(x) 



1 



1 



1 - X^l 1 - X^l 



=n 

a=l 



1 



1 - xy« 



A judicious introduction of additional variables allows us to actually compute quantities 
associated with the solutions of the diophantine equations (A.2 ). For example, if we want 
to get the solutions themselves, we can introduce a set of variables y^, oc G N, and take 



T'i(x;A/i,V2, 



1 



1 



1 — A^ix^i 1 — yix^^i 



n 



1 



The coefficient of x''^ is now a polynomial in the such that each monomial of the form 

= rii, Nvt = n2 . . . If we want to obtain 



tells us that there is a solution N. 



just the sum Y.a. l^kk' = ^ for all a G N. In this case, the coefficient of x''^ is 

a polynomial in the y such that the degree n of each monomial signals the existence of 
a solution to (A.2) with Y.a = ^- This idea has been used to obtain the generating 



function for the peaks given by (VI.6). Finally the substitution of the variables y^ by the 
specific numerical value yc = 2 would provide a generating function such that the coef- 
ficient of x'l^ would be ^sois(qi) n« 2^''", where the previous sum extends to the solutions 
of ( A.2[ > for a given qi and is zero if the equation has no solutions. 

As we really have to consider sets of decoupled diophantine equations the actual num- 
ber of solutions is the product of the number of solutions for each of them and hence the 
generating function becomes 



F(xi,x2,...)=n^^^^^)=nn 



1 



1 — xr 



i=1 i=l a=l 

Hence, for a = c\^^/p^ + c^iy/Vi + ■ ■ ■ / the number that we are looking for is 

Y_ 1 = [x?'xf ■■■]F(Xi,X2,...), 

cee{a) 

where [x^^ • ■ ■ ]F(xi , X2, . . .) denotes the coefficient of the x^^ x^"^ ■ • • term in the power 
series expansion of F(xi , X2, . . .) around (xi , X2, . . .) = (0, 0, . . .)• 

The factors appearing in the denominators of Dr(a) can be incorporated by modifying 
the previous generating function F^ as follows 



Hi(x) 



0! 



1! 



2! 



0! 



1! 



2! 



Now, it is immediate to obtain the number we are looking for as 

1 



[x^'xf • ■ ■]Fl(Xi,X2, . . .) , 



a=1 
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where 

oo oooo^ /oooo\ 

H(xi , . . .) = n = n n ^ =^^p\^y_ • (^.3) 

i=l i=l a=1 a=l / 

Finally, the factor in the numerator of Di-fa) is N!, with N = Nk^^- Iri order to 

introduce this factor we notice that the expansion of the exponential appearing in ( A.3 ) is 

(oooo \ / oo oo \ ^ / oo oo \ ^ 

^^xf" Ul+ ^^xf" +l ^^xf« +■■■ 
1=1 a=l / \ i=l a=l / ■ V i=l a=l / 

The key point now is to realize the N*-term in the above expansion. 



N! 



. i=l a=l 



gathers the contributions of those elements of C(a) satisfying precisely Y.i Y.a ^ 
This can be easily seen by i ntroducing the single variable y as above and checking that 
each term of the form ( A.4[ > appears now multiplied by y^. Then, it suffices to modify 
H(xi , X2, . . .) by multiplying each term in the previous expansion by N ! to get 



oooo \ /oooo \ /oooo 



G,(x,,x,,...) = l + ^^xf") + |^^xf'') +---+(LLxf«) +• 

\i=l a=l 

' oo oo 



1=1 a=l / \ i=l a=1 / \ i=l a=l 

oo oo 



i=l a=l 



Summarizing, given an area eigenvalue a = (]^^/p^ + c^iy/pi + ■ ■ • , its r-degeneracy is 
given by 

D,(a) = [x?'xf ■■■]G,(xi,X2,...). 



A.2. Generating functions for the m-degeneracy 



A straightforward extension of the argument given in the previous subsection to com- 
pute quantities associated with the solutions to ( A.2[ > by introducing auxiliary variables 
allows us to incorporate the projection constraint (or similar conditions) to explain the ra- 
degeneracy It is straightforward, in particular, to introduce the terms for the projection 
constraints (see subsections IIL3 IV.l| and IV.2[ > 




z — 



[ _1 ^. T oo oo / Vi- -1-1 — — 1 

[z — Z ) -I— r-i— r / Z —Z" 



Z — Z 



-1 



for DL , 
for GM , 
for ENP . 
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By doing this we get 

G°^(Z, Xi,X2, 
G™P(Z,X1,X2, 

V 1=1 ra=l 

and, from them, we obtain the black hole degeneracies as 

DC(a) = [z°][x?'x2^^---]GC(z,Xi,X2,...), 
where the index C refers to the counting scheme that we are interested in (i.e DL, GM or 




ENP). Notice that, as we mentioned in section VI the "master" generating functions G^ 
allow us to obtain other generating functions that can be used to study the behavior of 
the degeneracy for every conceivable subfamily of area eigenvalues. 



Appendix B: Explicit computations 

We give here a complete computation of the black hole degeneracy spectrum and en- 
tropy, corresponding to the first eigenvalues of the area spectrum, according to the DL 
prescription [21J. On one hand, we will provide explicit computations to concretely show 
how the different methods introduced in the main body of the paper work. On the other, 
we want to explore the behavior of the entropy for the smallest black holes. This will al- 
low us to see how the shapes of the black hole degeneracy spectrum -and of the entropy- 
as a function of the area arise. Important features such as the linear growth of the entropy 
with the area and the appearance of the periodicity observed in [9 , lOJ can be already seen 
at this level. It is also possible to directly study the role of the projection constraint in the 
definition of the entropy. The main steps of the computation are: 

1) Determination of all the area eigenvalues smaller than a fixed value a+. For practi- 
cal purposes we will consider areas smaller than 18 (in units of Anyl^). 

2) Computation of the black hole degeneracy spectrum by using generating functions. 

3) Determination of the entropy according to the Domagala-Lewandowski recipe. 



We will discuss these points with some detail. 

1) The density of area eigenvalues grows very quickly as a function of the area. There 
are 354 area eigenvalues smaller than 18. Though they can be easily handled by a com- 
puter, they are too many to be listed here so in the following we provide a table with only 
those eigenvalues smaller than 12. Nevertheless we will extend the plots that we give in 
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this appendix up to areas of around 18. For a given upper bound of the area a+, we have 
to find all the different numbers of the form 



1=1 



obtained by considering positive half -integers ji G N/2 and such that they are smaller 
than a+. In this process, different choices of the ji may actually give the same value. This 
source of degeneracy is taken into account in a precise way by our counting methods (the 
Pell equation, generating functions and so on). At this level we will just care about the 
area eigenvalues disregarding their degeneracies. 

For a+ = 18 the maximum allowed value of ji is obtained by solving the inequality 
A/j(j + 1) < 9; this gives j = 17/2. On the other hand, if each ji takes the smallest non-zero 



allowed value of 1/2, the maximum number of terms in the sum (B.l I is [18/V3J = 10. 
This means that we can generate all the sought for area eigenvalues by considering the 
points in the discrete set (N U {0}y^ contained in the ten-dimensional simplex defined by 
the condition < 17/2. The 1/10! factor reduction with respect to the computation 

extended to the cubic grid {1 /2, 1 , . . . , 1 7/2}^° is important to reduce the computing time. 

2) The generating functions that we have introduced in Appendix |A] depend, in prin- 
ciple, on an infinite number of variables Xi associated with the square-free numbers pi. 
For a finite subset of the area spectrum only the square roots of a finite number of them 
are relevant. This means that we will only have to consider the finite number of variables 
associated with them. If we take a < 1 8, the square roots of the squarefree numbers that 
appear in the area eigenvalues are listed in Table |lj The only variables that we need to 
write explicitly in the generating functions are Xi, X2, X3, X4, X5, X/, X9, Xio, Xig, X22, X27, Xgg, 
^119/ ^156/ ^198 • We have to consider then: 

G°^(Z, Xi , X2, X3, X4, X5, X7, X9, Xio, X18, X22, X27, Xgg, X119, Xi56, Xi98) 

= (1 -(z2 + z-2)xf -(zi6 + z-i^)x]2 -(z + z-^)x2 -(z^ + z-^)x^ 
-{z' + z-']xt -(z4 + z-4)x| -{z' + z-'Ws -{z' + z-']x^j 

-{z"' + Z-'']xt -(z3 + Z-3)Xio -(zl° + Z-l°)xfg -(z5 + Z-5)X22 

-{z'^ + Z-'VlJ -(Z"+Z-")X88 -(Zl3 + Z-I3)xii9 _(zl5 + ^-15)^^^^ _ ^^17 ^ ^-17)^^^^) 
G^J'^fXi , X2, X3, X4, Xs, X/, X9, Xio, Xi8, X22, X27, X88, Xi 19, Xisg, Xi98) 

= (1 -2xf -2x]^ -2x2 -2x^ -2x4 _2^2 _2^3 _2^3 

-2x^ -2xio -2x^g -2x22 -2x^7 -2Xg8 -2Xii9 -2Xi56 - 2Xi98)"V 

The coefficients [z°][x^^ • • ■x'^^l^]G^^ and [x^^ ■ ■ ■x1^l^]G^^, respectively, tell us the values 
of D°L(a) andD°L(a) for a = q^V2+- ■ •+qi98V323. For example, for a = 2^2+3 v^S+v^ 
(i.e. the row corresponding to n 71 in the table given at the end of this appendix) we 
get 



D^^[2V2 + 3V3 + yA5] = [z°][xfx^xio]G^^ = 120, 
D^^ilVl + 3V3 + v^) = [xfx^xio]G°^ = 640 . 



3) The entropy is obtained by adding up the values of D^^[a) up to a certain value of 
the area. 
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The table at the end of this appendix shows the explicit values of the black hole degen- 
eracies for the first area eigenvalues (corresponding to areas smaller than 12 in our units) 
and also their cumulative sum (that in the case of D^^[a) gives the exponential of the 
entropy). The results of an explicit computation up to areas of the order of 18 are shown 
in Figs. 13 to 15 



We want to make several comments at this point: First, in order to understand some 
of the features of the entropy and the role of the projection constraint, it is convenient to 
study the auxiliary entropy S< ( a) given by (III.27) and compare it with the actual entropy 
S< (a). Figure 13 shows both the exact values of S<(a) for all the area eigenvalues smaller 



than 1 8 and the values of S< (a) as a function defined for all areas a G [0, 1 8]. It is possible 
to see that, after a short transient regime, the entropy grows linearly with area and a 
characteristic staircase structure appears. As can be seen in the detailed plots of Fig. 14 
this is more evident for the largest areas considered, where the plot of the entropy can 
be effectively approximated by a smooth curve due to the increasing density of the area 
spectrum. The width of the steps is roughly 0.35 in our units. This is in satisfactory 
agreement with the prediction of equation ( VL4[ >. 



In order to obtain the physical entropy S< (a) the projection constraint must be incor- 
porated. It is instructive to compare the results obtained in this case with the ones de- 
scribed above. In particular Fig. [15] is the counterpart of Fig. 131 In this case the staircase 



structure is more evident because the width of the steps doubles (to a value around 0.7). 
This doubling is a consequence of the effective suppression of many configurations by 
the action of the projection constraint as justified in subsection VI.3 : these configurations 
are not shown in Fig. 15 Figure 15 shows also the values of the entropy corresponding 
to the prequantized values of the areas (that are just integer values in units of 47ty£p). For 
this subset of areas the growth of the entropy as a function of the area is linear (with the 
additional logarithmic corrections). The imprint of the staircase structure is this case is 
the presence of some larger than average jumps ("double jumps") in the value of the en- 
tropy for successive prequantized areas. Finally the last figure (Fig. [16]) shows the black 
hole degeneracy spectrum with its characteristic peak structure. As it can be seen the 
peaks are quite pronounced and the distance between them corresponds to the width of 
the steps in the entropy. 
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FIG. 14. Detail of Fig. 13 for areas a € [14, 18]. Notice the steps that appear for the area values 
considered in this plot. 




FIG. 15. Plot of the value of S< (a) for both the points in the area spectrum and all the values of 
the area smaller than 18 (in units of 47ty£p). We also show the entropy values corresponding to 
prequantized values of the area. 
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FIG. 16. Logarithmic plot of D^^[a) for all the points in the area spectrum smaller than 18 (com- 
pare this figure with Fig. [7]| and a detail of the last two peaks. The arrows mark the position of 
the bands as predicted by equation ( VI.4 1. Notice that there are many points in the area spectrum 
for which D^Lfal = 0. 
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Appendix C: A group theoretic treatment of the projection constraint 



As pointed out in subsection IV. 1 the problem of finding the m-degeneracy for the GM 



counting is equivalent to determining the number of irreducible representations -taking 
into account multiplicities- that appear in the tensor product 0^^^[ki/2]. The reason 
is that each of the irreducible representations that appear in the decomposition of the 
tensor product as a direct sum has, precisely, one basic state with zero total third spin 
component. 

In order to solve this problem we propose a solution based on techniques developed 
in the context of conformal field theories 11251 142| . The starting point is to write the tensor 
product of two SU(2) irreducible representations in the form 



ki 
2 



k2 
2 



k2 



k3=0 



k3 

2 



in terms of the fusion numbers ^k^k2- taking into account that the tensor product and 
direct sum of irreducible representations have a direct translation into the behavior of the 
characters, and in particular: 

1. The algebra of the characters of the SU(2) irreducible representations satisfies 



Xk, ■ Xk2 



ka 



•3 > 



2. and the characters of irreducible representations are ortonormal with respect to the 
SU(2)-scalar product induced by the (normalized) Haar measure, i.e. 



(XknXki) 



SU(2) - 



S3 



Xk,Xk2d|J.s3 = 6(ki,k2) 



we can easily obtain the number of irreducible representations in the composition of 
spin-k/2 associated with a configuration c = {(k, N^)} as 



d$;^(c) 



L(nxkSxk') 



SU(2) 



k'=0 k 



2 ~ s.n->llc+1)e^sm(K- + l)9 



sin 9 



(C.l) 



-Izl n 



7t n 

k'=0-^" 



sin^^ (k+ 1)9\ sin(k' + 1)9 



sin''^9 



sin 9 



sin^ 9d9 



Notice that the group coordinate 9 is naturally defined in 9 G [0, tt] but, due to the parity 
properties of the integrand in (C.l I, it is possible to extend the integration in 9 to [0, 2tc]. 
This will prove useful to solve the present problem in yet another different way. The 
key idea is to realize that we are actually composing SU(2) representations and asking 
how many states with vanishing total third spin component appear in such composition. 
This is done by computing, with the help of the characters of the representations, the 
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multiplicity of the 0-irreducible representation of U( 1 ) in the decomposition of the tensor 
product of the SU(2) representations involved. The characters rji^, k G Z, of the U(1) 
irreducible representations are orthonormal with respect to the standard scalar product 
in the circle 



(■fik, ,'ni<2)u(i) 



r2n 



S' 



rik,rik2 dfj-si 



g-ik,egik2e 



de 

Ik 



We can obtain the number that we are looking for just by projecting the product of char- 
acters of the SU(2) representations onto the character rji^^o of the U(1 ) irreducible repre- 
sentation 



qGMf 



Ufl) 



1 r27r 
271 



n 



/sin(k+ 1)9 
V sine 



d0. 



(C.2) 



The last expression coincides with the one derived in subsection IV.l by using generating 
functions and shows an interesting interplay between the counting of SU(2) and U(1) 
labels. 

Similar considerations allow us to obtain the formulas corresponding to the ENP and 
DL countings. In the first case we have to find the multiplicity of the singlet SU(2) ir- 
reducible representation in the composition of the representations appearing in a given 
configuration. This can be trivially obtained by projecting over the character Xo- In this 
way we get 



•271 



sm 



sin (k+ 1)6 
sin 9 



d9 

n 



Finally, the result for the DL counting can also be obtained by using reducible U(l ) repre- 
sentations with characters fjk = ilk + il-k and projecting over rjo. 



C(c) 



1 r27t 

In 



f3(2cosk9)^^de. 
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